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ABSTRACT 

This  work  focuses  on  the  construction  of  equilibrated  two-phase  antiplane 
shear  deformations  of  a  non-elliptic  isotropic  and  incompressible  hyperelastic  ma¬ 
terial.  It  is  shown  that  this  material  can  sustain  metastable  two-phase  equilibria 
which  are  neither  piecewise  homogeneous  nor  axisymmetric,  but,  rather,  involve 
non-planar  interfaces  which  completely  segregate  inhomogeneously  deformed  ma¬ 
terial  in  distinct  elliptic  phases.  These  results  are  obtained  by  studying  a  con¬ 
strained  boundary  value  problem  involving  an  interface  across  which  the  deforma¬ 
tion  gradient  jumps.  The  boundary  value  problem  is  recast  as  an  integral  equation 
and  conditions  on  the  interface  sufficient  to  guarantee  the  existence  of  a  solution 
to  this  equation  are  obtained.  The  contraints,  which  enforce  the  segregation  of 
material  in  the  two  elliptic  phases,  are  then  studied.  Sufficient  conditions  for  their 
satisfaction  are  also  secured.  These  involve  additional  restrictions  on  the  inter¬ 
face  across  which  the  deformation  gradient  jumps — which,  with  all  restrictions 
satisfied,  constitutes  a  phase  boundary.  An  uncountably  infinite  number  of  such 
phase  boundaries  are  shown  to  exist.  It  is  demonstrated  that,  for  each  of  these, 
there  exists  a  solution — unique  up  to  an  additive  constant — for  the  constrained 
boundary  value  problem.  As  an  illustration,  approximate  solutions  which  corre¬ 
spond  to  a  particular  class  of  phase  boundaries  are  then  constructed.  Finally,  the 
kinetics  and  stability  of  an  arbitrary  element  within  this  class  of  phase  boundaries 
are  analyzed  in  the  context  of  a  quasistatic  motion. 
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1.  INTRODUCTION 

Finite  elastic  equilibria  with  discontinuous  deformation  gradients  have  fig¬ 
ured  prominently  in  recent  continuum  mechanical  treatments  of  displacive  solid- 
solid  phase  transformations.  Models  of  this  sort  are  pertinent  to  the  investiga¬ 
tion  of  shape  memory,  twinning  and  transformation  toughening  in  solids.  Shape 
memory  and  transformation  toughening  occur  in  metallic  and  ceramic  alloys, 
respectively,  while  twinning  can  be  found  in  both  metallic  and  ceramic  alloys. 
Micrographs  of  multiphase  equilibrium  states  in  alloys,  such  as  those  presented 
by  Zackay,  Justusson,  &  Schmatz  [27]  and  Porter  &  Heuer  [21],  often  dis¬ 
play  configurations  wherein  the  various  phases  are  segregated  by  geometrically 
complicated  interfaces.  One  question  which  arises  regarding  the  aformentioned 
continuum  mechanical  idealizations  of  such  materials  is  whether  they  axe  capable 
of  capturing  the  morphological  complexity  of  such  deformations.  As  a  first  step 
toward  answering  this  question,  this  work  focuses,  within  the  context  of  a  par¬ 
ticular  class  of  hypothetical  materials,  on  the  construction  of  equilibria  involving 
coexistent  phases  segregated  by  surfaces  which — although  not  as  morphologically 
complex  as  those  displayed  in  [21,27] — are,  at  least,  non-planar. 

In  a  homogeneous,  hyperelastic  material  discontinuous  deformation  gradients 
occur  only  if  the  relevant  elastic  potential  allows  for  a  loss  of  ellipticity,  at  certain 
values  of  the  deformation  gradient,  in  the  associated  displacement  equations  of 
equilibrium.1  Materials  characterized  by  elastic  potentials  which  allow  such  a 
loss  of  ellipticity  are  referred  to  as  non-elliptic.  Of  particular  importance  in  this 
work  are  non-elliptic  materials  which  have  at  least  two  disjoint  elliptic  phases. 
Examples  of  such  materials  are  provided  by  Ericksen  [12]  in  the  context  of  a 
one- dimensional  bar  theory,  and  by  Abeyaratne  [3]  in  his  study  involving  a  spe¬ 
cial  class  of  incompressible,  isotropic  materials.  Abeyaratne  [3],  Ball  &  James 
[10],  Gurtin  [16],  and  Silling  [25]  have  demonstrated  that  materials  of  this  sort 
support  equilibrium  states  which  display  coexistent  elliptic  phases  and,  in  addi- 

1  For  a  discussion  of  this  issue  see,  for  instance,  Rosakis  [24], 
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tion,  minimize  the  relevant  energy  functional.  As  a  result  of  the  latter  property 
these  states  are  referred  to  as  mechanically  stable.  The  associated  deformation 
fields  in  all  of  the  foregoing  works  are  either  piecewise  homogeneous  or  axisym- 
metric.  In  the  case  of  equilibrated  piecewise  homogeneous  deformations  the  as¬ 
sociated  phase  boundaries  must  be  planar.  Ball  &  James  [10]  and  Silling  [25] 
have  shown,  however,  that  energy-minimizing  sequences  of  piecewise  homogeneous 
mechanically  stable  deformations  may  possess  limits  which  are  metastable  as  op¬ 
posed  to  mechanically  stable  and,  moreover,  involve  non-planax  phase  boundaries. 
On  the  other  hand,  ABEYARATNE  [2]  and  SlLLING  [25]  have  constructed,  respec¬ 
tively,  asymptotic  and  numerical  solutions  to  a  boundary  value  problem  involving 
a  mode  III  crack  in  a  particular  subclass  of  incompressible,  isotropic  non-elliptic 
materials.  These  solutions  are  not  mechanically  stable  and  axe  neit  er  piece 
wise  homogeneous  nor  axisymmetric;  furthermore,  they  include  the  non-elliptic 
material  phase  and,  in  addition,  transitions  between  the  elliptic  and  non-elliptic 
phases  which  do  not  involve  jumps  in  the  deformation  gradient.  These  solutions 
do,  however,  involve  surfaces  which  separate  the  two  elliptic  phases  present  in  the 
deformation.  The  relevant  interfaces  are,  moreover,  non-planar.  Rosakis  [23]  has 
recently  shown  that  a  special  class  of  anisotropic  non-elliptic  materials  is  capable 
of  sustaining  equilibria  in  which  a  family  of  cusped  lenticular  inclusions  of  one 
elliptic  phase  reside  in  a  matrix  of  another  elliptic  phase.  These  states  are,  in 
general,  metastable. 

As  yet  there  are  no  results  pertaining  to  the  existence,  in  non-elliptic  isotropic 
hyperelastic  materials,  of  multiphase  equilibrium  states  which  are  free  of  the  non- 
elliptic  phase  and  are  neither  piecewise  homogeneous  nor  axisymmetric.  The 
primary  objective  of  this  investigation  is  to  prove  constructively  that  a  class  of 
non-elliptic  isotropic  incompressible  hyperelastic  materials  is  capable  of  sustaining 
deformations  of  this  type.  These  deformations  will  typically  be  metastable — like 
those  associated  with  the  limits  of  the  aformentioned  minimizing  sequences  of 
piecewise  homogeneous  deformations  and  the  states  constructed  by  Rosakis  [23]. 
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Although  non-planar  phase  interfaces  may,  in  reality,  reflect  anisotropic  effects, 
these  results  show  that  they  can  exist  within  the  context  of  a  model  which  does 
not  take  anisotropy  into  consideration.  Isotropic  materials  may,  consequently,  be 
useful  in  preliminary  studies  of  the  kinetics  and  stabilty  of  interfaces  between 
phases.  These  issues  are  taken  up  briefly  in  the  final  section  of  this  work  and, 
more  thoroughly,  by  FRIED  [13]  in  a  linear  stability  analysis  of  states  involving 
planar  phase  interfaces  for  a  class  of  non-elliptic  isotropic  materials. 

Chapter  2  is  devoted  to  preliminaries.  After  a  brief  overview  of  the  notation 
to  be  used,  Section  2.1  introduces  the  kinematics  and  fundamental  balance  prin¬ 
ciples  which  will  be  needed  in  the  following.  Section  2.2  explains  the  constitutive 
restrictions  which  will  be  adhered  to  throughout  this  work.  Section  2.3  begins  by 
introducing  the  concept  of  a  quasistatic  motion.  It  then  discusses  the  notions  of 
mechanical  dissipation  and  driving  traction  associated  with  surfaces  across  which 
the  deformation  gradient  jumps;  these  lead  naturally  to  the  consideration  of  a 
kinetic  relation  and  the  associated  kinetic  res-ponse  function.  In  the  final  section 
of  Chapter  2,  the  kinematics  are  specialized  to  those  of  antiplane  shear. 

Chapter  3  focuses  upon  the  solution  of  a  particular  constrained  boundary 
value  problem,  in  antiplane  shear,  involving  the  field  equations  and  jump  con- 
di'ions  put  forth  in  Section  2.4.  After  formulating  the  problem  in  Section  3.1, 
a  representation  for  the  solution  of  the  boundary  value  problem  is  presented  in 
Section  3.2.  This  representation  is  indeterminate  in  that  it  involves  the  unknown 
jump  in  the  normal  derivative  of  the  displacement  field  over  an  interface  across 
which  the  deformation  gradient  is  discontinuous.  In  Section  3.3  an  integral  equa¬ 
tion  is  derived  for  the  unknown  jump  in  the  normal  derivative  of  displacement 
in  terms  of  a  parameterization  of  the  interface.  Sufficient  conditions  for  the  ex¬ 
istence  of  a  unique  solution  of  this  integral  equation  are  then  obtained.  These 
constitute  analytical  restrictions  on  the  interface  geometry.  It  transpires  that 
there  exist  an  uncountably  infinite  number  of  interfaces  which  comply  with  these 
resrictions.  It  is  then  shown  that  for  each  of  these  interfaces  there  exists  a  sclu 


tion,  unique  up  to  an  additive  constant,  tc  the  boundary  value  problem  stated 
in  Section  3.1.  In  Section  3.4  the  constraints  which  enforce  segregation  of  the 
elliptic  phases  are  analyzed.  These  impose  further  analytical  restrictions  on  the 
interface  geometry.  It  is  shown  that,  within  the  set  of  interfaces  which  allow  a 
solution  to  the  boundary  value  problem,  there  exist  an  uncountably  infinite  num¬ 
ber  of  interfaces  which  also  satisfy  these  restrictions  and,  hence,  allow  a  solution 
to  the  constrained  boundary  value  problem.  Each  of  these  solutions  involves  a 
non-planar  and  non-axisymmetric  phase  interface  which  separates  elliptic  phases 
subjected  to  inhomogeneous  deformations.  Section  3.5  illustrates  the  results  of  the 
two  preceeding  sections  in  determining  a  particular  class  of  non-planar  surfaces 
for  which  the  constrained  boundary  value  problem  can  be  solved.  Approxima¬ 
tions  for  the  strain  and  displacement  fields  corresponding  to  the  solutions  of  the 
appropriate  family  of  constrained  boundary  value  problems  are  then  constructed. 

The  last  chapter  is  concerned  with  the  kinetics  and  stability  of  slowly  propa¬ 
gating  phase  boundaries.  In  Section  4.1  the  distribution  of  driving  traction  along 
a  phase  interface  of  the  kind  constructed  in  Section  3.5  is  calculated.  Section 
4.2  is  concerned  with  observations  pertaining  to  the  kinetics  of  such  a  surface. 
Ingredients  crucial  to  this  analysis  are  the  kinetic  relation  and  response  function 
introduced  in  Section  2.3.  This  section  concludes  with  results  that  relate  the 
monotonicity  of  the  kinetic  response  function  for  a  particular  material  and  the 
kinetic  stability  of  that  material.  These  final  results  are  consistent  with  those 
obtained  in  [13]. 
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2.  PRELIMINARIES 

2.1.  Notation,  kinematics  and  balance  principles.  In  the  following  IR 
and  (T  denote  the  sets  of  real  and  complex  numbers.  The  intervals  (0,oo)  and 
[0,oo)  are  represented  by  1R+  and  1R+.  The  symbol  IRn,  with  n  equal  to  2  or  3, 
represents  real  n-dimensional  space  equiped  with  the  standard  Euclidean  norm. 

-  O 

If  U  is  a  set,  then  its  closure,  interior  and  boundary  are  designated  by  U ,  U,  and 
dU ,  respectively.  The  complement  of  a  set  V  with  respect  to  U  is  written  as 
U  \  V.  Given  a  function  ip  :  U  — >  W  and  a  subset  V  of  U ,  ip(V)  stands  for  the 
image  of  V  under  the  map  ip. 

Vectors  and  linear  transformations  from  Ft3  to  1R3  (referred  to  herein  as 
tensors)  are  distinguished  from  scalars  with  the  aid  of  boldface  type — lower  and 
upper  case  for  vectors  and  tensors,  respectively.  The  symbol  £  refers  to  the  set 
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of  tensors,  £+  denotes  the  set  of  all  tensors  with  positive  determinant,  and  S 
stands  for  the  collection  of  all  symmetric  positive  definite  tensors.  The  set  of  unit 
vectors  in  1R3  is  designated  by  M.  If  F  is  in  £  then  FT  represents  its  transpose; 
if,  moreover,  det  F  ^  0,  then  the  inverse  of  F  and  its  transpose  are  written  as 
F-1  a.id  F-T,  respectively.  When  component  notation  is  used,  Greek  indices 
range  only  over  {1,2};  summation  of  repeated  indices  over  the  appropriate  range 
is  implicit. 

Let  q  lie  in  [1,  oo).  Then  a  function  ip  :  1R  —*  JR  is  an  element  of  Lq{JR)  if  it 
is  q  integrable  on  IR — that  is,  if  its  Lq  norm  over  IR, 

+  oo 

=  (  J  MONO*. 

—  OO 

is  defined.  Similarly,  ip  is  an  element  of  L°°(1R)  if  ip  is  bounded  on  IR — that  is,  if 
its  norm  over  IR, 

=  SUP  1^(01. 

teiR 

exists.  A  function  :  IR 2  — *  IR  is  an  element  of  Lq(IR2)  or  L°°(IR2),  respectively, 
if  the  analogous  Lq  or  L°°  norm  over  IR2  exists. 
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Consider,  now,  a  body  B  which,  in  a  reference  configuration,  occupies  a  region 
TZ  contained  in  Ft3.  Let  the  invertible  mapping  y  :  TZ  — * ►  TZm,  with 

y(x)  =  x  +  u(x)  Vx  6  TZ,  (2.1.1) 

characterize  a  deformation  of  B  from  the  reference  configuration  onto  a  config¬ 
uration  that  occupies  the  region  7 Z»  in  HI3.  Assume  that  the  deformation  y,  or 
equivalently  the  displacement  u,  is  continuous  and  possesses  piecewise  continuous 
first  and  second  gradients  on  7Z.  Let  5  be  the  set  of  points  contained  in  TZ  defined 
so  that  y  is  differentiable  on  the  set  TZ\S.  Introduce  the  deformation  gradient 
tensor  F  :  R  \  S  — >  £+  by 


F(x)  =  Vy(x)  Vxe7 Z\S,  (2.1.2) 

and  assume  that  the  associated  Jacohian  determinant,  J  :  TZ  \  S  — ►  ffl,  of  y  is 
strictly  positive  on  its  domain  of  definition: 

J(x)  =  detF(x)  >  0  Vx<=7Z\S.  (2.1.3) 

+ 

The  left  Cauchy- Green  tensor  G  :  TZ  \  S  — >  S  corresponding  to  the  deformation 
y  is  given  by 

G(x)  =  F(x)Fr(x)  Vx  €  TZ  \  S.  (2.1.4) 

The  deformation  invariants  associated  with  y  exist  on  TZ  \  S  and  are  supplied 
through  the  fundamental  scalar  invariants  of  G: 

ii(G)  =  tr  G,  72(G)  =  |((trG)2-tr(G2)),  /3(G)  =  detG.  (2.1.5) 

With  the  above  kinematic  antecedents  in  place  introduce  the  nominal  mass 
density  p  :  TZ  — *  HZ+,  the  nominal  body  force  per  unit  mass  b  :  TZ  — *  M3,  and  the 


nominal  stress  tensor  S  :  7v  \  S  — >  C,  and  suppose  that  p  and  b  are  continuous 
on  7Z,  while  S  is  continuous  on  1Z  \  S,  piecewise  continuous  on  7 and  also  has  a 
piecewise  continuous  gradient  on  7 Z.  In  the  absence  of  constitutive  assumptions 
relating  the  stress  to  the  deformation  gradient,  the  sets  over  which  S  and  F  suffer 
jumps  need  not  be  equivalent.  The  scope  of  this  investigation  is  limited,  however, 
to  elastic  materials  wherein  stress  is  continuously  related  to  strain — hence.  S,  like 
F,  is  continuous  on  7 Z\S.  Let  p,  be  the  mass  density  in  the  deformed  configuration 
associated  with  y.  Given  a  regular  subregion  V  of  7 Z,  let  m  :  dV  — ►  A”  denote 
the  unit  outward  normal  to  dV .  Then  the  global  balance  laws  of  mass,  and — in 
the  absence  of  inertia — force  and  moment  equilibrium  require  that 


JpdV=  j  p.  d\\ 

V  y(T’) 

J  Sm  dA  +  J  ph  dV  =  0, 

dV  V 

j  y  A  Sm  dA  +  J  y  A  pb  dV  =  0, 


(2.1.6) 


(2.1.7) 


(2.1.8) 


respectively,  for  every  regular  subregion  V  contained  in  71. 

Localization  of  the  balance  laws  (2.1.6)— (2.1.8)  at  an  arbitrary  point  con¬ 
tained  in  the  interior  of  7 Z  \  S  yields  the  following  familiar  field  equations: 


P  =  on  7 Z\S, 

V  ■  S  +  ph  =  0  on  7Z  \  S, 
SFr=FSr  on  7v  \  5, 


(2.1.9) 


Suppose,  from  now  on,  that  the  set  5  is  a  regular  surface.  Then,  localization 
of  (2.1 .6)— (2. 1 .8)  at  an  arbitrary  point  in  5  yields  the  following  jump  conditions: 


!My)Jl  =  0  on  s, 

[Sn]  =  0  on  5, 


(2.1.10) 
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where,  given  a  generic  field  quantity  g  :  1Z  \  S  — ►  1R  which  jumps  across  5,  [<?]  is 
defined  through 

[<7(x)J  =  lim  (fl>(x  +  /m(x))  -  g(x  -  hn(x)))  Vx  6  S,  (2.1.11) 

h  \,0 

with  n  :  5  — ►  AT  a  unit  normal  to  5.  Observe  from  the  jump  condition  (2. 1 . 1 0 ) i 
that  the  mass  denity  in  the  deformed  state  p„  is  only  defined  on  y(77.  \  S).  Equa¬ 
tions  (2.1.9)i  and  (2.1. 10) i  are,  evidently,  completely  decoupled  from  equations 
(2.1. 9)2,3  and  (2.1.10)2;  that  is,  given  a  solution  to,  say,  a  boundary  value  prob¬ 
lem  involving  the  latter  set  of  equations,  p„  can  be  calculated  a  posteriori.  For 
this  reason  equations  (2. 1 .9)  1  and  (2.1.10)]  will  be  disregarded  in  the  subsequent 
analysis.  In  addition  to  the  jump  conditions  given  in  (2.1.10),  the  stipulated 
continuity  of  y  gives 

[u]  =  0  on  5.  (2.1.12) 

2.2.  Constitutive  assumptions.  Let  B  be  composed  of  a  hyperelastic  ma¬ 
terial  which  is  homogeneous,  isotropic  and  incompressible.  Since  B  is  hyperelastic 
its  mechanical  response  is  governed  by  an  elastic  potential  or  strain  energy  per 
unit  reference  volume.  The  homogeneity  of  B  implies  that  the  elastic  potential 
does  not  depend  explicitly  on  position  in  the  reference  configuration.  Further¬ 
more,  because  B  is  isotropic  the  elastic  potential  can  depend  on  the  deformation 
gradient  F  only  through  the  deformation  invariants  Ik{G )  defined  in  (2.1.7).  The 
incompressibility  of  B  requires  that  the  deformation  y  be  isochoric ,  *.  e. , 

73(G(x))  =  J2(x)  =  1  Vxg7 l\S.  (2.2.1) 

An  additional  consequence  of  isotropy  is,  therefore,  that  the  elastic  potential  can 
be  expressed  as  a  function  solely  of  the  first  two  deformation  invariants.  It  can 
also  be  demonstrated  via  (2.1.5)  that,  when  (2.1.1)  holds,  /Q(G(x))  >  3  for  all 
x  contained  in  7 Z\S.  Now,  let  W  :  [3,  00)  x  [3,  00)  — *  1R  denote  an  elastic 
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potentiai  which  characterizes  B  and  assume  that  W  is  continuously  differentiable 
with  piecewise  continuous  second  derivatives  on  its  domain  of  definition.  The 
nominal  stress  response  of  B  is  then  determined  through  W  up  to  an  arbitrary 
pressure  p  :  Tl  \  S  —>  Ft  required  to  accomodate  the  kinematic  constraint  (2.2.1) 
imposed  by  the  incompressibility  of  B:  viz., 

S  =  2(W/1(/)F  +  W/j(/)(/1(G)1-G)g)  -PF~t  on  7 l\S,  (2.2.2) 

where  1 :  71  \  S  — ►  [3,  oo)  x  [3,  oo)  is  given  by 

/(x)=  (;,(G(x)),/2(G(x)))  Vx€R\S. 

Following  Gurtin  [16],  let  the  class  of  generalized  neo-Hookean  materials  re¬ 
fer  to  that  subset  of  hyperelastic  materials,  first  introduced  by  Knowles  [19], 
which  are  homogeneous,  isotropic  and  incompressible  with  elastic  potential  in¬ 
dependent  of  the  second  deformation  invariant  (2. 1.5)2.  Assume,  henceforth, 
that  B  is  composed  of  a  generalized  neo-Hookean  material  with  elastic  potential 
W  :  [3,  oo)  — ►  JR,  where  W  is  continuously  differentiable  with  piecewise  contin¬ 
uous  derivative  on  [3,  oo).  Then,  by  (2.2.2),  the  nominal  stress  response  of  B  is 
determined  by 


S  =  2lK'(/1(G))F-pF-r  on  H\S.  (2.2.3) 

Suppose  also  that  the  elastic  potential  is  normalized  so  that 

W(3)  =  0.  (2.2.4) 

Choose  a  rectangular  Cartesian  frame  X  =  {0;ei,e2,e3}  and  consider  the 
response  of  the  material  at  hand  to  a  simple  shear  deformation  y  given  by 


y(x)  =  x  +  7(x  •  e2)e!  Vx  £  U, 


(2.2.5) 
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where  the  constant  7 — assumed  non-negative  without  loss  of  generality — denotes 
the  amount  of  shear.  From  (2.1.2),  (2.2.3)  and  (2.2.5)  the  nominal  shear  stress 
corresponding  to  the  deformation  y  is,  for  each  7  in  JR+,  found  to  be 

ei  •  Se2  =  27W'(3  4-  ■y2)  =:  t(j).  (2.2.6) 

In  [19-20]  Knowles  demonstrates  that  the  components  of  nominal  and 
Cauchy  shear  stress  are,  in  the  present  setting,  equal.  The  function  r  :  1R+  — >  M 
is,  hence,  referred  to  as  the  shear  stress  response  function  of  the  generalized 
neo-Hookean  material,  characterized  by  W,  in  simple  shear.  An  immediate  con¬ 
sequence  of  (2.2.4)  and  (2.2.6)  is 


V  b  — 3 

w(h)=  j  VIi  6  [3,oo),  (2.2.7) 

0 

so  that  the  response  of  a  generalized  neo-Hookean  material,  in  all  three  dimen¬ 
sional  deformations,  is,  up  to  a  hydrostatic  pressure,  completely  characterized 
by  specifying  the  shear  stress  response  function  r.  Define  the  secant  modulus  in 
shear  M  :  1R+  —*  JR  of  a  generalized  neo-Hookean  material  with  elastic  potential 
W  by 

M(j)  =  2W'(3  +  72)  V7  e  M+,  (2.2.8) 

and  assume  that,  in  compliance  with  the  Baker- Ericks en  inequality, 

M(7)  >0  V7  €  R+.  (2.2.9) 

Assume,  also,  that  M( 0)  >  0  so  that  the  infinitesimal  shear  modulus  of  the 
material  at  hand  is  positive.  Note  from  (2.2.6)  and  (2.2.8)  that  the  shear  stress 
response  function  r  must  also  satisfy 


r(0)  =  0,  r'(0)  =  M(  0). 


(2.2.10) 
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Observe,  also,  that  the  stipulated  smoothness  of  W  guarantees  that  both  r  and 
M  are  piecewise  continuously  differentiable  on  1R+. 

It  is  worth  remarking  that,  despite  the  significant  restrictions  which  have  been 
placed  upon  the  class  of  materials  which  will  be  considered  in  this  investigation, 
no  presuppositions  have  been  made  regarding  the  sign  of  the  derivative — where  it 
exists — of  the  shear  stress  response  function  corresponding  to  the  generalized  neo- 
Hookean  material  through  (2.2.6).  In  [20]  Knowles  shows  that  the  monotonicity 
of  the  shear  stress  response  function  r  is  related  directly  to  the  ellipticity  of  the 
generalized  neo-Hookean  material  which  it  characterizes:  if  r  is  not  a  monotone 
increasing  function  on  its  domain  of  definition  then  the  associated  material  is  non- 
elliptic.  This  investigation  will  make  use  of  a  particular  subclass  of  non-elliptic 
generalized  neo-Hookean  materials,  first  suggested  by  Abeyaratne  [3];  this  class 
of  materials  is  characterized  by  the  set  of  shear  stress  response  functions  r  which 

are  continuous  on  2R+  and  piecewise  continuously  differentiable  on  1R.+  \  {7,7}, 

* 

where  0  <  7  <  7,  such  that 

* 

r' >  0  on  jR+\[7,7],  r' <  0  on  (7,7).  (2.2.11) 

*  * 

The  sets  of  shear  strains  lying  in  the  intervals  [0,7)  and  (7, 00)  are  referred  to  as 

* 

the  high  and  low  strain  phases  of  the  generalized  neo-Hookean  material  specified 
by  the  shear  stress  response  function  r.  Together  the  high  and  low  strain  phases  of 
such  a  material  comprise  its  elliptic  phases.  A  generalized  neo-Hookean  material 
characterized  by  a  shear  stress  response  function  of  this  type  will  be  referred 
to  herein  as  a  three-phase  material.  See  Figure  1  for  a  graph  of  a  shear  stress 
response  function  typical  of  those  which  specify  three-phase  materials.  Within 
the  class  of  three-phase  materials  special  attention  will  be  given  those  materials 
(proposed  by  Abeyaratne  in  [2])  for  which 

M17  if  7  e  [0,7], 

* 

d( 7)  if  7  €[7,7], 

* 

M27  if  7  £  [%°°), 


tp(  7)  = 


(2.2.12) 
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where  the  function  d  :  [7,7]  — >  1R  is  linear  in  its  argument.  Observe  that,  in 

* 

accordance  with  (2.2.11)2,  d  is  required  to  decrease  on  (7,7).  A  further  conse- 

* 

quence  of  (2.2.11)  is  that  /zi  must  be  greater  than  [i 2  which  must  itself  be  positive. 
Figure  2  shows  the  graph  of  rp. 

2.3.  Dissipation,  driving  traction  and  the  kinetic  relation.  For  the 

purposes  of  this  section  it  is  necessary  to  consider  a  one  parameter  family  of 
deformations  y(-,t)  :  %  — *  TZt  where  t,  which  denotes  time,  increases  from  to  to 
1 1.  It  is  assumed  that  y(x,-)  is  continuous  with  piecewise  continuous  first  and 
second  derivatives  for  each  fixed  x  in  %.  Let  St  be  a  regular  surface,  with  unit 
normal  n(-,f)  :  S<  — ►  AT,  contained  in  7 Zt  for  each  value  of  t  in  [to,fi].  The  fields 
u  (■,/)  :  7£  -  IR3,  F(-,f)  :K\St-+  £+,  b(-,f)  :  TZ  — ►  JR3,  and  S(-,t)  :  7 Z\St^C 
are,  at  each  f  contained  in  [fo,ti],  the  obvious  counterparts  of  those  introduced  in 
Section  2.1.  A  one  parameter  family  of  deformations  of  this  sort  is  referred  to  as 
a  quasistatic  motion  if  the  above  quantities  and  the  nominal  mass  density  satisfy 
the  field  equations 


V  •  S(-,f)  4- pb(-,f)  =  0  on  7 Z\St  Vf  G  [<o,*i], 

S(-,()Fr(  ,()  =  F(-,()ST(-,0  on  fc\St  V<  €[<„,(!), 


(2.3.1) 


the  jump  condition 


[S(-,On(-,01  =  0  on  St  Vt€[«o,<  1],  (2.3.3) 

and  the  kinematic  condition  of  displacement  continuity 

|u(-,f)l  =  0  on  St  Wt€[t0,ti).  (2.3.4) 

Knowles  [18]  has  shown  that,  in  a  quasistatic  motion,  the  presence  of  a 
moving  surface  of  discontinuity  St  of  the  type  considered  here  has  an  effect  on 
the  balance  of  mechanical  energy.  Let  V  be  a  regular  subregion  contained  in  1Z. 


L 
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In  [18]  it  is  demonstrated  that  the  difference  in  the  rate  of  work  of  the  mechanical 
forces  external  to  V  and  the  rate  at  which  energy  is  stored  in  V  is  given  by 

Sa(t-V)=  !  f(x,t)Vn(x,t)dA  Vie  Mi],  (2.3.5) 

stnv 

where,  for  each  t  in  [fo,fi],  /(•,<)  :  St  — +  IR  is  the  scalar  driving  traction  and 
Vn(’,t )  :  St  — *  IR  is  the  normal  velocity  of  the  interface  (in  the  reference  configu¬ 
ration).  The  function  £,(•;  V)  :  [fo,<i]  — ♦  IR  is  referred  to  as  the  rate  of  dissipation 
of  mechanical  energy  associated  with  the  region  V.  It  has  been  shown  by  Yatomi 
&  Nishimura  [26]  as  well  as  Abeyaratne  &  Knowles  [7]  that  the  form  of  the 
driving  traction  for  a  hyperelastic  material  is,  in  the  quasistatic  setting,  supplied 

by 

/(■,()  =  [^(F(-,())]-S(-,()  [F(-,i)]  on  S,  Vie  [(„,(,],  (2.3.6) 

+  — 

where  S(-,t)  (reap.,  S(-,t))  is  the  limiting  value  of  the  field  S(-,i)  on  the  side  of 

the  interface  into  which  the  unit  normal  n (•,<)  is  (reap.,  is  not)  directed  at  t  in 

Mi]. 

When  treated  from  a  thermomechanical  perspective,  the  dissipation  rate  can 
be  shown  to  be  identical  to  the  product  of  the  temperature  and  the  rate  of  entropy 
production — provided  that  the  temperature  is  spatially  uniform  and  independent 
of  time.2  The  Clausius- Duhem  inequality  then  requires  that  the  dissipation  rate 
associated  with  a  quasistatic  motion  of  the  kind  envisioned  here  be  non-negative, 
that  is 

0  Vfe[fo,*i],  (2.3.7) 

for  every  regular  subregion  V  contained  in  TZ.  A  localization  of  (2.3.5)  at  an 
arbitrary  point  on  the  interface  therefore  yields  the  inequality 

f{-,t)Vn(',t)  >  0  on  St  Vie  Mi]  (2.3.8) 


2 


For  a  detailed  discussion  of  the  these  issues  see  Abeyaratne  &  Knowles  [7], 
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as  a  condition  imposed  for  the  admissibility  of  the  quasistatic  motion. 

In  the  context  of  a  motion  which  involves  such  an  interface  it  is  necessary  (see 
[1]  and  [3-8])  to  supplement,  in  some  fashion,  the  constitutive  information  which 
relates  the  stress  and  strain  fields.  An  approach  to  this  pioneered  by  Abeyaratne 
&  Knowles  [7]  entails  the  provision  of  a  kinetic  relation  which  gives  the  normal 
velocity  of  the  interface  in  terms  of  the  driving  traction  which  acts  thereon  or  vice 
versa.  In  the  former  case  one  specifies  a  function  V  :  1R  — +  JR  such  that 

vn  =  V(f)  V/  6  M.  (2.3.9) 

Here  V  is  referred  to  as  the  kinetic  response  function.  If  the  function  V  is  such 
that  V(f)f  >  0  on  JR  then  (2.3.8)  is  automatically  satisfied  and  the  kinetic  re¬ 
sponse  function  is  itself  referred  to  as  admissible.  If  an  admissible  kinetic  response 
function  is  continuous  on  M,  then  it  must  satisfy  V’(O)  =  0.  If,  in  addition  to  be¬ 
ing  admissible,  V  is  continuously  differentiable  on  JR,  then  V'(0)  >  0.  Otherwise 
admissibility  implies  nothing  with  regard  to  the  sign  of  the  derivative  of  a  smooth 
kinetic  response  function  V .  All  kinetic  response  functions  considered  herein  are 
assumed  to  be  admissible. 

In  the  work  of  Abeyaratne  [3],  Ball  &  James  [10],  Gurtin  [16],  Gurtin  & 
Temam  [15],  and  Silling  [25]  the  necessary  additional  constitutive  information  is 
provided  by  setting  the  driving  traction  equal  to  zero  on  St  for  all  t  in  [to  >  ^ l ]  •  This 
amounts  to  prescribing  a  particular  mte  independent  kinetic  relation  whereby 
energy  is  conserved;  it  is,  furthermore,  a  necessary  consequence  of  requiring  that 
a  suitable  energy  functional  be  minimized  at  each  t  in  [to ,  ti]  (see  Abeyaratne 
[!])• 

2.4.  Antiplane  shear  of  a  generalized  neo-Hookean  material.  Sup¬ 
pose,  from  now  on,  that  It  is  a  cylindrical  region  and  choose  a  rectangular  Carte¬ 
sian  frame  X  =  {0;ei,e2,e3}  so  that  the  unit  base  vector  e3  is  parallel  to  the 
generatrix  of  It.  The  deformation  y  defined  through  (2.1.1)  consists  of  an  an- 
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tiplane  shear  if  it  is  of  the  form 

y(x)  =  x-f  u(ii,x2)e3  Vx  €  (2.4.1) 

That  is,  the  displacement  field  intrinsic  to  an  antiplane  shear  deformation  has 
only  one  nonzero  component  in  the  e3  direction  which  is  independent  of  the  x3- 
coordinate.  In  (2.4.1)  xQ  =  x  •  e0  for  each  x  contained  in  TZ.  The  function  u  will 
be  referred  to  as  the  out-of-plane  displacement  field.  Inspection  of  (2.4.1)  reveals 
that  any  discontinuities  in  the  gradient  of  y  must  be  due  to  discontinuities  in  the 
out-of-plane  displacement  field  and,  hence,  occur  across  surfaces  which  do  not 
vary  with  the  x3-coordinate. 

Knowles  [19]  has  demonstrated  that,  although  not  every  hyperelastic, 
isotropic  and  incompressible  material  can  sustain  antiplane  shear  deformations, 
all  generalized  neo-Hookean  materials  are  capable  of  doing  so.  It  has  been  shown 
(Knowles  [19-20])  that  for  such  materials  the  local  balance  equations  (2. 1.9)2, 3 
reduce,  in  the  absense  of  body  forces  and  under  the  asssumption  that  the  nominal 
stress  tensor  is  independent  of  the  X3-coordinate,  to  the  scalar  equation 

(Af(7)tt,a),«  =  0  on  V\C.  (2.4.2) 

Here  M  is  the  secant  modulus  in  shear  as  defined  in  (2.2.8),  7  :  V\C  — >  1R  is  the 
shear  strain  field  given  by 

7(xi,x2)=  Vu>a(xi,X2)u,a(xi,x2)  V(ii, x2)  €  P  \  C,  (2-4.3) 

V  is  a  plane  region  with  shape  determined  by  a  generic  cross  section  of  H,  and 
C  is  a  curve  contained  in  T>  and  determined  similarly  by  a  cross  section  of  the 
surface  across  which  the  deformation  gradient  jumps.  Furthermore,  the  jump 
condition  (2.1.10)2  reduces,  for  a  generalized  neo-Hookean  material  subjected  to 
antiplane  shear,  to 


[M(j)u,a  n0J  =  0  on  C, 


(2.4.4) 
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where  n  :  C  — >  Af  is  a  unit  normal  to  C,  while  (2.1.15)  becomes 

[uj  =0  on  C.  (2.4.5) 

It  is  also  readily  shown  that  the  driving  traction  /,  introduced  in  Section  2.3, 
for  a  generalized  neo-Hookean  material  subjected  to  an  antiplane  shear  deforma¬ 
tion  involving  a  discontinuity  in  the  gradient  of  displacement  across  a  curve  C  is 
given  by 

/  =  [W(3  +  72)]  —  M(j)u,a  |u,a  ]  on  C.  (2.4.6) 

In  (2.4.6)  u,a  and  u,a  refer  to  the  limiting  values  of  the  gradient  of  the  out-of- 
plane  displacement  field  on  the  side  of  the  curve  C  into  which  and  out  of  which 
the  unit  normal  n  points,  respectively.  Evidently,  7  and  7  are  given  in  terms  of 
u,a  and  u,Q  by  (2.4.3). 
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3.  STUDY  OF  A  CONSTRAINED  BOUNDARY  VALUE  PROBLEM 
IN  THE  ANTIPLANE  SHEAR  OF  A  THREE-PHASE  MATERIAL 

This  chapter  focuses,  in  the  context  of  antiplane  shear,  on  the  construction  of 
two-phase  equilibria  of  a  body  composed  of  the  particular  non-elliptic  generalized 
neo-Hookean  material  with  shear  stress  response  function  rp  defined  via  (2.2.12). 
These  equilibria  will  involve  non-planar  interfaces  which  segregate  material  in 
different  elliptic  phases.  The  interfaces  will  be  described  by  surfaces  Qs  of  the 
form 

Qa  =  {x  6  JR3  |  x\  =  s(x2),x 2  G  1R,  x 36  IR} 

where  s  is  twice  continuously  differentiable  on  IR,  s ^  is  in  L°°(IR)  fl  L2(IR)  for 
n  —  0, 1,2,  and 

lim  s(x  2)  =  0. 

i2-±oo  v  ’ 

In  Section  3.1  a  boundary  value  problem  for  the  out-of-plane  displacement  field 
associated  with  two-phase  antiplane  shear  deformations  of  a  three-phase  material 
is  formulated  and  specialized  to  the  case  of  the  material  with  shear  stress  response 
rp.  This  boundary  value  problem  is  supplemented  by  a  set  of  constraints  which 
require  that  the  non-elliptic  phase  of  the  relevant  material  is  absent  and,  moreover, 
that  the  elliptic  phases  are  segregated.  In  Section  3.2  the  boundary  value  problem 
is  converted  into  an  integral  equation  for  the  jump  in  the  normal  derivative  of 
the  out-of-plane  displacement  field  across  Q9-  In  Section  3.3  it  is  shown  that 
there  exists  a  unique  solution  to  this  integral  equation  for  every  Qs  defined  by  a 
function  s  which,  in  addition  to  the  restrictions  delineated  above,  satisfies 


+00 

(  ( |«'(*2)||."(*,)|dr2)* 

J  ’  Ml  —  1*2 

—00 


where  /ij  and  fi2  are  the  moduli  associated  with  the  elliptic  phases  of  the  material 
defined  by  rp.  It  is  then  demonstrated  that  for  each  such  s  there  exists  a  unique 
(up  to  an  arbitrary  additive  constant)  solution  to  the  boundary  value  problem 
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stated  in  Section  3.1.  Since,  however,  the  constraints  mentioned  above  are  not 
necessarily  satisfied  by  any  of  these  solutions,  the  deformation  associated  with 
a  given  solution  does  not  necessarily  constitute  a  two-phase  equilibrium  state  of 
the  type  sought  after  here.  In  Section  3.4  it  is  shown  that,  provided  a  certain 
functional  of  s  is  sufficiently  small — in  a  sense  to  be  made  precise — then  there 
exists  a  unique  solution  to  the  constrained  boundary  value  problem  and,  hence,  a 
two-phase  equilibrium  state  of  the  type  sought  after  here.  The  concluding  section 
of  this  chapter  is  concerned  with  the  construction  of  a  class  of  two-phase  states 
which  involve  non-planar  interfaces  separating  material  in  distinct  elliptic  phases. 

3.1.  Formulation  and  reduction  of  the  boundary  value  problem 
and  phase  segregation  requirements.  Suppose  that  B  is  composed  of  a 
three-phase  material  and  that  the  cylinder  71  degenerates  to  occupy  all  of  ]R3. 
Let  the  rectangular  Cartesian  frame  X  be  as  in  Section  2.4.  Consider  the  effect 
of  subjecting  71  to  a  particular  antiplane  loading  whereby,  independent  of  the 
^-coordinate,  the  shear  strain  approaches  uniform  values  of  7/  as  ij  tends  to 
—00  and  7r  as  X\  tends  to  +00.  Assume  that  7 /  is  greater  than  7  and  that 

7r  lies  strictly  between  0  and  7;  the  prescribed  remote  shear  strains  associated 

* 

with  the  loading  are,  thus,  in  the  high  and  low  strain  phases  of  the  material  at 
hand  as  x\  approaches  —00  and  +00,  respectively.  If  7/  and  yr  are  chosen  so 
that  the  corresponding  remote  shear  stresses  7(7/ )  and  r(7r)  are  equal  then — for 
every  three-phase  material — there  exists,  modulo  an  arbitrary  additive  constant,  a 
unique  one  parameter  family  of  pairwise  homogeneous  out-of-plane  displacement 
fields  ua :  1R  —>  ]R  which  satisfy  the  equilibrium  equation  in  (2.4.2)  on  JR 2  \  Ca, 
with  the  straight  line  Ca  given  by  {xQe0  E  IR2 |xi  =  a,x2  6  JR},  the  jump 
conditions  in  (2.4.4)  and  (2.4.5)  on  Ca  and,  of  course,  the  decay  requirements 
associated  with  the  prescribed  conditions  at  x\  —  ±00.  The  function  ua  is  given 
by 


,  v  f  7/(xi  “  a )  if  xi  <  a, 

0  1  (  7r(xi  — a )  if  >  a. 


(3.1.1) 
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Here  a  determines  the  point  of  intersection  of  the  plane  surface  Qa,  given  by 
Ca  x  2R,  with  the  xi-axis.  Note  that,  for  each  fixed  a,  the  pairwise  homogeneous 
deformation  associated  with  (3.1.1)  through  (2.4.1)  involves  exclusively  the  elliptic 
phases  of  the  material  under  consideration  and  that  these  are  segregated  by  Qa; 
the  deformation  associated  with  (3.1.1)  will  be  referred  to  as  a  globally  elliptic 
pairwise  homogeneous  equilibrium  state.  The  interface  Qa  associated  with  such  a 
state  will,  in  turn,  be  referred  to  as  a  phase  boundary.  Observe  that  the  qualitative 
character  of  the  equilibrium  state  associated  with  ua  is  clearly  unaffected  by  the 
value  of  a. 

Envision  a  generalization  of  the  globally  elliptic  pairwise  homogeneous  equi¬ 
librium  state  wherein  the  kinematics  remain  those  of  antiplane  shear  and  the 
loading  conditions  are  as  described  at  the  outset  of  this  section  but  the  planar 
phase  boundary  is  replaced  by  a  non-planar  interface  Qs  with  cross  section  C9 
where,  for  simplicity, 

Qs=C,xR  (3.1.2) 

with 

C„  =  {zQeQ  €  R2\xi  =  s(x2),X2  G  R}.  (3.1.3) 

Assume  that  the  state  is  equilibrated  in  the  sense  that  the  balance  equation  in 
(2.4.2)  holds  on  R2  \  Cs  while  the  jump  conditions  in  (2.4.4)  and  (2.4.5)  are 
satisfied  on  Cs.  Clearly,  if  such  a  state  exists,  the  deformation  field  intrinsic  to  it 
must  be  inhomogeneous  on  either  side  of  the  interface  Qs.  Observe  that  even  if  a 
three-phase  material  is  capable  of  sustaining  a  deformation  of  this  kind  the  shear 
strain  field  may  not,  in  general,  be  distributed  so  that  only  the  elliptic  phases  of 
the  material  are  present;  if,  however,  this  is  the  case  and,  furthermore,  the  high 
and  low  strain  phases  of  the  relevant  material  are  segregated  by  the  interface  Qs 
then  the  deformation  will  be  said  to  constitute  a  globally  elliptic  inhomogeneous 
two -phase  equilibrium  state  with  phase  boundary  Qs. 

Consider,  now,  the  geometry  of  the  curve  C,  which  determines  the  phase 
boundary  Qt  essential  to  a  globally  elliptic  inhomogeneous  two-phase  equilibrium 
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state.  Since  the  shear  strain  field  is  constant  as  x\  approaches  ±oo,  it  is  clear 
that  s  must  be  bounded  on  JR  in  order  for  Ca  to  qualify  as  a  cross  section  of  the 
phase  boundary  Qa.  The  kinematics  and  boundary  conditions  place  no  further 
restrictions  on  the  geometry  of  Ca . 

Observe  that  if,  in  addition  to  being  bounded  and  continuous  on  1R,  s  satisfies 
one  or  both  of 

lim  s(r2)  =  c,  lim  5(2-2)  =  c,  (3.1.4) 

z  2— ♦  — 00  rj- >+00 

where  c  and  c  are  real  constants,  then  the  loading  must  be  restricted  so  that  the 
far  field  shear  stresses  r(7/)  and  r( 7r)  are  equal.  To  see  this  suppose  that  (3.1  -4)i 
holds.  Then,  as  12  approaches  —00  the  phase  boundary  becomes  planar  and  the 
local  character  of  the  deformation  begins  to  resemble  a  pairwise  homogeneous 
state.  Since  the  far  field  shear  strains  7/  and  7r  Eire  constant  the  local  shear 
strains  must  match  these  appropriately  on  either  side  of  Ca  as  x\  approaches  —00. 
Hence,  the  local  shear  stresses  must  match  their  far  field  counterparts  r(7/)  Emd 
r( 7r)  and,  by  the  jump  condition  in  (2.4.4)  which  holds  on  Ca,  r(-yi)  =  r( 7,.). 
A  completely  analogous  argument  can  be  constructed  if  (3.1.4)2  holds  instead  of 
(3.1  -4)i .  Certainly,  if  both  of  (3.1.4)  hold,  the  result  is  still  true.  Note,  however, 
that  if  neither  of  (3.1.4)  hold,  and,  hence,  the  curve  C,  is  merely  bounded,  there  is 
no  reason  to  rule  out — a  priori — loading  conditions  wherein  the  far  field  stresses 
are  unequal. 

Assume,  henceforth,  that  (3.1.4)  holds  with  c  and  c  equal  to,  say,  c.  Recalling 
the  rol*’  of  a  in  (3.1.1),  there  is  certainly  no  additioned  loss  in  generality  incurred 
by  talcing  c  =  0.  In  this  case  (3.1.4)  becomes 

lim  s(x  2)  =  0. 

ij  — *±oo 

Let  U  be  the  set  of  functions  defined  by 

U  =  {t/>  :  JR  — ►  JR  |  ip  €  C(JR),  lim  ^(*2)  =  0,  t/>  ^  0  on  JR}. 

r  j— »±oo 


(3.1.5) 
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Assume,  henceforth,  that  s  is  an  element  of  the  set 

A  =  u  nvnw, 


(3.1.6) 


where  V  and  W  are  defined  by 

V  =  {</>  :  R  -  R\xl>  G  C2(Jff?),V’(n)  G  L2(lR),n  =  0,1,2},  (3.1.7) 

and 

W={xP.n->  JR  10  €  C2(J?),^(n)  G  L°°(]R), n  =  0,1,2},  (3.1.8) 


respectively. 

Given  an  element  s  of  A  which  describes  an  interface  Q3  it  is  convenient  to 
define  plane  sets  T>[  and  VTS  by 

V\  =  {xQea  €  1R2\xx  <s(x2),*2  €^},  £>;  =  JR2  \  T>[.  (3.1.9) 


Cleaxly,  the  union  and  intersection  of  Vl3  and  V3  form  generic  cross-sections 
of  the  cylinder  71  and  the  phase  boundary  Qs,  respectively.  Note,  also,  that 
if  s  is  an  element  of  A  then,  by  (3.1.3)  and  its  assumed  smoothness,  a  unit 
normal  to  Q,  exists  everywhere  on  Qs  and  depends  only  on  the  X2-coordinate. 
Let  n  :  IR  — »  M  designate  the  unit  normal  to  Qs  which  points  into  the  region  of 

O 

low  strain — VI  x  IR.  Then  the  representation  for  n  is  computed  easily  from  the 
definitions  of  Qj  and  C,  and  is  given  by 


n(r2) 


ei  -  s'(x2)e2 
yj\  +  s'(x2)2 


Vx2  G  IR. 


(3.1.10) 


Now,  if  a  three-phase  material  is  capable  of  sustaining  a  globally  elliptic 
inhomogeneous  two- phase  equilibrium  state  of  antiplane  shear  with  phase  bound¬ 
ary  Qt  then  the  out-of-plane  displacement  field  u  associated  with  the  deformation 
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through  (2.4.1)  must,  by  virtue  of  (2.4.2),  (2.4.4)  and  (2.4.5),  satisfy  the  following 
field  equation  and  jump  conditions: 

(M(y)u,a),a  =  0  on  JR2  \CS, 

lM(y)u,QnQJ  =  0  on  C3,  (3.1.11) 

fu]  =  0  on  C3, 

with  n  as  indicated  in  (3.1.10);  in  order  to  comply  with  the  prescribed  loading  it 
suffices  to  require  that  the  gradient  of  u  satisfies  the  following  asymptotic  decay 
conditions: 


^io  (•*■  1  1  ')®<* 


f  7/ei  +  o(l)  as  xi  — >  —  oo, 
l7rei+o(l)  as  xj  — ►  +oo, 


(3.1.12) 


moreover,  in  order  to  assure  that  only  the  elliptic  phases  of  the  material  at  hand 
are  present  and  are  segregated  by  Qs,  the  shear  strain  field  7,  given  in  terms  of 
the  gradient  of  u  by  (2.4.3),  must  conform  to  the  following  inequalities: 


7e(7>°°)  on  £»>  7  £  [0, 7)  on  (3.1.13) 

* 

where  Vls  and  V?  are  given  by  (3.1.9).  These  inequalities  will  be  referred  to  as 
the  phase  segregation  requirement. 

Given  a  three-phase  material,  (3.1.11)— (3.1.12)  comprise,  for  each  fixed  s 
contained  in  A ,  a  boundary  value  problem  in  the  out-of-plane  displacement  field 
u,  while  (3.1.13)  acts  as  a  system  of  constraints  thereon.  Together  (3.1.11)— 
(3.1.13)  will  be  referred  to  as  the  constrained  boundary  value  problem  in  u  for 
the  three-phase  material  with  secant  modulus  in  shear  M.  Given  a  particular 
three-phase  material  the  constrained  boundary  value  problem  need  not  have  a 
solution  for  any  function  s  in  A.  The  study  of  (3.1.11)— (3.1.13)  for  a  specific 
material  may,  however,  serve  as  a  means  to  determine  a  subset  of  A  for  which 
the  constrained  boundary  value  problem  is  soluble. 
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Before  proceeding  note  that  the  jump  conditions  (3.1.11)2,3  holding  across 
C.  can  be  recast  to  read 


M(tW-)+.-))|(»(-)+,-)  =  M(rW-)-.-))|(*(-)-,  •)  on  2R, 
u(s(-)+,-)  =  u(s(-)-,-)  on  IR, 


where  the  +  and  —  symbols  indicate  the  limiting  values  of  the  appropriate  quan¬ 
tities  on  the  high  and  low  strain  sides  of  the  interface,  respectively. 

For  simplicity  attention  will,  for  the  remainder  of  this  work,  be  restricted 
to  the  constrained  boundary  value  problem  for  the  material  characterized  by  the 
shear  stress  response  function  rp  defined  in  (2.2.12).  In  this  case  the  form  of 
the  shear  stress  response  function  is  such  that  the  secant  modulus  in  shear  M 
is  constant  in  both  the  high  and  low  strain  elliptic  phases;  hence,  (3.1.1 1  )i  and 
(3.1.14)i  reduce  to 


u,aQ  =  0  on  IR 2  \  CS , 


(3.1.15) 


The  analytical  difficulties  of  the  special  constrained  boundary  value  problem  posed 
by  (3.1.15),  (3.1.14)2,  (3.1.12)  and  (3.1.13)  are  certainly  less  daunting  than  those 
encountered  in  the  analogous  problem  for  a  more  general  three-phase  material. 
For  each  fixed  s  in  A  the  only  non-linearity  which  encumbers  the  problem  asso¬ 
ciated  with  tp  is  that  imposed  by  the  strain  constraints  (3.1.13).  In  the  present 
absence  of  results  pertaining  to  the  existence  of  globally  elliptic  two-phase  equi¬ 
libria  in  arbitrary  three-phase  materials,  any  results  which  can  be  obtained  for 
this  particular  material  constitute  progress  toward  a  qualitative  understanding  of 
the  more  general  issue. 

As  a  first  step  in  analyzing  the  constrained  boundary  value  problem  com¬ 
prised  by  (3.1.15),  (3.1.14)2,  (3.1.12)  and  (3.1.13)  it  is  convenient  to  introduce  a 
reduced  out-of-plane  displacement  field  v  :  ]R2  \  CB  — >  IR  specified  via 


v(xux2)  =  u(xux2)  -  u0(xux2)  V(x!,x2)  €  lR2\Ct, 


(3.1.16) 
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where  u0  :  IR2  \  Ca  — ►  IR  is  furnished  by 


Uq(x\,X2)  =  H  (s(x2)  -  +  H  (x  1  -  s(x2  ))  ~irX\  V(xi,X2)  G  2R2  \C3, 


and  H  :  1R  \  {0}  — ►  iR  is  the  Heavyside  function: 


0  if 
1  if 


xi  <  0, 

xi  >  0. 


Solving  for  u  in  (3.1.16)  and  inserting  the  result  in  (3.1.15),  (3.1.14)2  and 
(3.1.12)  shows,  with  the  aid  of  the  definitions  of  u0  and  H,  that  the  reduced  out- 
of-plane  displacement  field  v  must  satisfy  the  following  boundary  value  problem: 


V,aa  =  0  on  lR2\Ca, 

=  /*2fri^ °n 

u(s(-)+»-)  -»(«(•)“*•)  =(7/-7r>  on  iR, 

w,a  (xi, -)eQ  =  o(l)  as  xx  — ♦  ±oo  on  iR. 


(3.1.17) 


Note  that  in  deriving  (3.1.17)2  use  has  been  made  of  the  equality  of  remote  shear 
stresses — which,  as  shown  at  the  beginning  of  this  section,  is  a  necessary  conse¬ 
quence  of  (3.1.5).  The  phase  segregation  requirement  (3.1.13)  can  be  written — 
after  appropriate  substitution  for  u — in  terms  of  the  components  of  the  gradient 
of  v  as  follows: 

7  ^  V,a  Via  +7/(2u,i  +7/)  on  Z>J, 

(3.1.18) 

0  <  V,a  V,a  +Jr(2v,i  +Jr)  <  7 2  On  Vf . 

* 

The  boundary  value  problem  in  (3.1.17)  will  be  referred  to  as  the  reduced 
boundary  value  problem  with  the  implicit  understanding  that  it  is  in  the  reduced 
out-of-plane  displacement  field  v  and  for  the  special  three-phase  material  with 
shear  stress  response  function  rp.  The  system  of  inequalities  in  (3.1.18)  will  be 
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labelled  the  reduced  phase  segregation  requirement.  It  is  clear  from  the  simple 
relation  between  the  primitive  out-of-plane  displacement  field  u  and  its  reduced 
counterpart  v  that  any  solution  to  the  reduced  problem  yields  a  solution  to  the 
original  problem.  The  next  section  will  focus  on  obtaining  a  representation  for 
the  solution  to  the  reduced  boundary  value  problem  with  the  reduced  phase  segre¬ 
gation  requirement  held  in  abeyance.  This  representation  will  lead  to  an  integral 
equation  which,  for  each  fixed  s  in  A ,  may  be  analyzed  in  place  of  the  associated 
reduced  boundary  value  problem. 

3.2.  Reformulation  of  the  reduced  boundary  value  problem  as  an 
integral  equation.  Let  s  be  an  arbitrary  element  of  A.  Since,  by  (3.1.17)i, 
the  reduced  out-of-plane  displacement  field  v  is  harmonic  on  M2  \  C9  the  jump 
conditions  (3.1.17)2  and  (3.1.17)3  suggest  that  v  can  be  represented,  modulo  an 
arbitrary  additive  constant,  as  the  sum  of  a  single-  and  a  double-layer  potential 
along  the  curve  C3.3  The  densities  of  the  appropriate  single-  and  double- layer 
potentials  are  given,  respectively,  in  terms  of  the  jumps  in  the  normal  derivative 
of  v  and  of  v  itself  across  the  curve  C„.  From  the  jump  condition  (3.1.17)3  it  is 
clear  that  the  density  of  the  double-layer  potential  must  be  given  by  (7 /  —  7r)s  on 
1R.  Since,  by  the  definition  of  the  shear  stress  response  function  rp,  the  moduli 
and  fj.2  which  appear  in  (3.1.17)2  are  required  to  be  unequal,  this  jump  condition 
does  not  yield  direct  information  regarding  the  form  of  the  density  of  the  single¬ 
layer  potential.  It  is,  therefore,  necessary  to  designate  the  jump  in  the  normal 
derivative  of  v  across  C$  in  terms  of  an  unknown  function — say  (7/  —  7 r)(f>,  where 
it  is  assumed,  until  demonstrated  otherwise,  that  <f>  :  1R  — ►  1R  does  not  vanish 
identically  on  1R.  Hence,  the  proposed  representation  for  the  reduced  out-of-plane 
displacement  field  v  takes  the  form 

v{xux2)  =  7<  7r  (5^(ri,r2)  +  x2))  V(xj,x2)  G  2R2 \C,,  (3.2.1) 

Ztt 

where  the  functions  :  JR2  — »  1R  and  Dt  :  1R2  \Ca  — ►  2R  issue,  respectively, 

3  For  an  overview  of  the  relevant  potential  theory  see  CoURANT  &  HILBERT  [11]. 
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from  the  single-  and  double-layer  potentials  on  Ca  with  densities  <j>  and  s  and  are 
given  by 


+oo 

S<t>(x  l,X2)  =  J G{(xi,X2,€)4>(£)y/l  +  s'(£)2  d£  V(xi,x2 )  €  1R2,  (3.2.2) 

—  OO 


D,(x  i,x2)=  Jg32(Xi,x2,Z  )s(t)d(  V(xi  ,x2)elR2\Cs.  (3.2.3) 

—  OO 

The  kernels  G{  :  ( 2R 2  \  CS)  x  M  1R  and  Gj  :  1R2  x  1R  — ►  M  which  appear  in 
(3.2.2)  and  (3.2.3)  are  given,  for  each  £  contained  in  1R,  by 


<??(*!, *2,0  =  In  y/(x!  -  s(0)2  +  (X2  ~  £)2  V(*i,x2)  £  lR2\Ca,  (3.2.4) 


and 


G‘2(xux2,0 


(Sl  ~  3(0)  ~  (X2  ~  £)$'(£) 

(xi  -  «(£))2  +  (x2  -  02 


v(xl5x2)  e  JR2. 


(3.2.5) 


Consider,  now,  the  issue  of  verifying  the  status  of  the  representation  (3.2.1)  as 
a  solution  to  the  boundary  value  problem  (3.1.17).  Since  the  single-  and  double- 
layer  potentials  are  harmonic,  by  construction,  on  1R2  \  Ca,  it  is  evident  that  the 
function  v  given  by  (3.2.1)-(3.2.5)  satisfies  (3.1.17)i.  A  series  of  direct  calcu¬ 
lations  too  long  to  display  here  show  that  v,i  (xi,-)  and  v,2  (*j,-)  both  behave 
asymptotically  like  0(l/xi)  as  xi  approaches  ±00  on  1R  so  that  the  represen¬ 
tation  (3.2.1)-(3.2.5)  complies  with  (3.1.17)4  and,  hence,  the  loading  conditions. 
Since  the  single-layer  term  (3.2.2)  is  continuous  on  1R2  and  the  double-layer  term 
(3.2.3)  has  been  constructed  so  that  is  possesses  a  jump  of  2 its  across  the  curve 
Ca  it  is  also  clear  that  (3.2.1)-(3.2.5)  furnishes  a  representation  of  v  which  sat¬ 
isfies  the  jump  condition  in  (3.1.17)3.  The  only  remaining  requirement  which 
must  be  satisfied  by  (3.2.1)-(3.2.5)  in  order  for  it  to  provide  a  solution  to  the 
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reduced  boundary  value  problem  is  the  jump  condition  (3.1.17)2  involving  the 
normal  derivative  of  v.  A  straightforward  but  tedious  calculation  using  standard 
results  from  potential  theory  delivers  the  limits  of  the  normal  derivative  of  v  on 
either  side  of  Ca  in  the  form 


+  oo 

t; Ww)±'i,) = - 

—  oo 
+  oo 

— OO 


±  ^~^<t>{x2)  VX2  € 


(3.2.6) 


where,  for  each  fixed  Z2  in  1R,  I»(x2,-)  :  1R  \  {2:2}  — ►  1R  and  Ks(x2,  ■):  M  —*  JR 
are  given,  repectively,  by 


I»(X2,0 


(${x2)  -  3(Q)a,(g2)  +  (J2  ~  Q 

(s(i2)-5(0)2  +  (x2  -02 


V£€2R\{*2}, 


(3.2.7) 


and 


K(r  ^  (a(*2)-«(0)-(*2-0<r(*2)  r  m 


(3.2.8) 


Observe  from  (3.2.8)  that  as  £  approaches  Ia(x2,  £)  is  singular  for  each 
a?2  in  JR  in  that 


J»(I2,£)~ - 7  as  £  -»  *2  Vi2  £  iR. 

x2  -  £ 

Hence,  the  integral  involving  IB  in  (3.2.6)  must,  as  indicated,  be  taken  in  the  sense 
of  the  Cauchy  •principal  value.  One  also  finds  that,  as  x2  approaches  £,  Is(x2,£) 
is  singular  for  each  £  in  1R  in  a  manner  entirely  analogous  to  that  displayed 
above.  On  the  other  hand,  an  examination  of  (3.2.6)  reveals  that  the  behavior  of 
K„(x2,£)  as  either  £  approaches  a  fixed  x2  in  JR  or  as  x2  approaches  a  fixed  £  in 
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2R  is  regular;  in  fact,  since  s  is  an  element  of  A ,  the  limits 


lim  Ka(x2,£)  = 
5— *12 


*“(x2) 

2(l  +  s'(x2)2) 


Vr2  G  2R, 


and 


lim  Ka(x2, 0  = 


*w(0 

2(1  +s'(02) 


Aem, 


both  exist  and  are  finite. 

Recall,  now,  that  the  function  <f>  which  appears  in  the  second  term  on  the  right 
hand  side  of  (3.2.6)  is  unknown.  The  jump  condition  (3.1.17)2  serves,  therefore,  as 
a  device  by  which  this  function  can  be  determined.  An  appropriate  substitution  of 
(3.2.6)  into  (3.1.17)2  yields — after  collecting  terms  and  dropping  a  non- vanishing 
common  factor — the  following  equation: 


(fXl  +  /i2)<£  + 


A*i  ~  1^2 
■K^/l  +  (s')2 


+  oo 


— oo 


+  A02dZ 


~  H2 

TT-y/l  -f  (s')2 


+oo 


on 


2R. 


(3.2.9) 


Observe  that  (3.2.9)  constitutes,  for  each  fixed  s  contained  in  A,  a  linear  integral 
equation  to  be  solved  for  <f>  on  2R.  The  integral  equation  in  (3.2.9)  can  be  simplified 
by  making  a  few  modest  substitutions;  toward  this  end  define  <p  :  1R  — ►  M  by 


<p(x2)  =  <t>(x2)y/l  +  s'(x2)2  Vx2  G  2R,  (3.2.10) 


and  introduce  a  real  constant  A  through  the  relation 

(3.2. 

Recall  from  the  definition  of  rp  that  the  moduli  //j  and  fi2  satisfy  0  <  \i2  <  /zj ;  A 
must,  consequently,  lie  strictly  between  0  and  1/tt.  Continuing  with  the  simplifi¬ 
cation  of  (3.2.9),  multiply  and  then  divide  both  sides  of  the  integral  equation  by 


^1  ~  ^2 
7r(^l  +  W) 
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yj\  +  (s')2  and  (^i  +  ^2)1  respectively,  to  obtain,  with  the  aid  of  the  definitions 
(3.2.10)  and  (3.2.11)  the  following  alternative  to  (3.2.9): 

+00  +00 

+  =  on  R.  (3.2.12) 

—  OO  —  OO 

For  the  purpose  of  facilitating  the  forthcoming  discussion  introduce,  for  each 
function  s  contained  in  A,  an  operator  M a  such  that,  for  each  function  ip  the 
function  M.sip  is  given  by 

+00 

M»iP=  fK,(-,t)tP(t)dt  on  JR.  (3.2.13) 

—  OO 

In  addition,  let  a  function  f9  :  M  — >  JR  be  defined  for  each  s  in  A  via 

+00 

/,=  on  R.  (3.2.14) 

—  OO 

With  the  aid  of  (3.2.13)  and  (3.2.14),  (3.2.12)  can  be  recast  to  read 

<p  +  \  M$<p  =  on  R.  (3.2.15) 

Evidently  a  solution  <p  to  (3.2.15)  provides,  through  (3.2.10),  a  solution  to 
(3.2.9).  However,  it  is  also  clear  from  (3.2.1)-(3.2.3)  and  (3.2.10)  that,  given  (p, 
v  can  be  obtained  directly  in  the  form 

+00 

v(xi,i2)=  7r  J{G\{x l,X2,CM0  +  G|(Xl,X2,05(0)^ 

— OO 

V(*i,x2)  €  R2\Ca,  (3.2.16) 

which  obviates  the  need  to  consider  <p.  Hence,  for  each  s  in  A  the  task  of 
constructing  a  solution  v  to  the  corresponding  reduced  boundary  value  problem 
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(3.1.17)  is  altered,  via  potential  theory,  to  one  of  constructing  a  solution  to  the 
corresponding  integral  equation  (3,2.15).  The  task  of  the  next  section  is  to  deter¬ 
mine  a  set  of  conditions  upon  s — in  addition  to  requiring  that  it  to  be  an  element 
of  A — which  are  sufficient  to  guarantee  the  existence  of  a  solution  to  (3.2.15). 

3.3.  Analysis  of  the  integral  equation.  Suppose  that  s  is  contained 
in  A  and  consider  the  kernel  Ka  associated  with  it  by  (3.2.8).  Observe  that 
the  stipulated  smoothness  of  s  implies  that  Ka  is  a  continuous  function  on  JR  . 
Moreover,  since 

l*.(«„OI  V(*2,f)  (3.3.1) 

(Z2~C) 

the  boundedness  of  s"  on  JR  and  Taylor’s  theorem  imply  the  following  global 
estimate  for  the  modulus  of  Ka: 

|K.(**,OI<5«upU"(OI  =  5U"Il-(«)  v(x2,0eK2.  (3.3.2) 

2  £€jR  2 

Hence,  the  kernel  Ka  corresponding  to  any  s  in  A  is  continuous  and  bounded  on 
JR;  furthermore,  the  bound  is  given  explicitly  in  terms  of  a  functional  of  s — the 
L°°  norm  of  s"  over  1R.  If  the  integral  equation  held  over  a  compact  domain 
then  the  bound  (3.3.2)  would  lead,  for  each  fixed  A  in  (0, 1  /tt),  to  sufficient  con¬ 
ditions  in  terms  of  the  size  of  ||s"|[ioc(;R)  which  would  allow  the  construction  of 
a  unique  solution  to  the  integral  equation  via  a  uniformly  convergent  Neumann 
series.  Since  the  integral  equation  in  (3.2.15)  holds  over  JR,  it  will  be  convenient 
to  determine  conditions  on  functionals  of  s  other  than  its  L°°  norm  which  are  suf¬ 
ficient  to  guarantee  an  analogous  result.  Toward  this  end  consider  the  Neumann 
series  for  this  integral  equation.  This  series  is  readily  obtained  via  the  method  of 
successive  substitutions  and  is  given  by  #  :  1R  — ►  JR  as  defined  below: 

OO 

#  =  A  £(~A)"A<?/,  on  JR. 

n=0 


(3.3.3) 
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Observe  that,  with  the  aid  of  (3.2.13)  and  a  formal  interchange  of  summation  and 
integration,  #  satisfies 


#  +  \MS$  =  \fa  on  ]R. 

That  is,  provided  the  formal  operations  performed  above  can  be  justified,  $  fur¬ 
nishes  a  solution  to  the  integral  equation.  If  the  Neumann  series  converges  uni¬ 
formly  then  this  is  certainly  the  case.  Consider  the  following  geometric  series: 

OO 

9  —  ^  •  (3.3.4) 

n=0 

If  Ka  and  fa  are  elements  of  L2(1R2)  and  L2(Ft),  respectively,  and  the  L 2  norm 
of  Ka  over  Ft2  satisfies 

^  1-^"*  lz,2(iR2)  <  1  (3.3.5) 

then  (3.3.4)  will  converge.  Note  that  the  Neumann  series  is  majorized  by  the 
geometric  series.  Conditions  sufficient  to  guarantee  the  convergence  of  (3.3.4) 
are,  accordingly,  sufficient  to  assure  that  the  Neumann  series  converges  uniformly 
on  its  domain  of  definition  and,  therefore,  as  alluded  to  above,  that  <P  supplies  a 
solution  to  (3.2.15).  Provided  these  sufficient  conditions  are  in  force,  the  operator 
Ma  is,  moreover,  a  Fredholm  integral  operator  with  domain  and  range  L2(Ft). 
Hence,  the  Fredholm  alternative  holds  and  it  can  be  shown  that  the  solution  <£  to 
the  integral  equation  provided  by  the  Neumann  series  is  unique.4 

At  present  the  aforementioned  sufficient  conditions  are  only  of  value  if  there 
exist  functions  s  in  the  set  A  defined  by  (3.1.6)-(3.1.8)  for  which  they  hold.  It 
will  now  be  demonstrated  that  the  first  two  conditions  are  satisfied  for  every  s 
in  A  and  that  the  third  holds  for  every  s  contained  in  the  proper  subset  J  of  A 
defined  by 

+0O 

l  =  {seA\\yJ%(  J\s'(x2)\\s"(x2)\dx2)1'  <1}.  (3.3.6) 

—  oo 


4 


See  GARABEDIAN  [14]  for  a  discussion  of  Neumann  series  and  the  foregoing  results. 
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First  suppose  that  s  is  an  element  of  A.  Show  that  the  kernel  Ka  must 
consequent ially  be  square  integrable  on  its  domain  of  definition.  Let  ka  :  Ml2  — >  Ml 
be  defined  by 


(3.3.7) 


Note,  from  (3.3.1),  (3.3.7)  that 


+  oo  +oo 


+  00  +oo 


//  K]{x2,£)didx 2  <  II 


|s(x2)  -  s{i)  -  (X2  -  Qs'(x2)\2 


—  oo  — oo 


(*2  -  04 

—  oo  —oo 
-l-oo  -f  oo 

=  J  J k2s(x2,v)drjdx2. 


d(  dx  2 


(3.3.8) 


Hence,  to  demonstrate  that  Ka  is  contained  in  L2(Ml2)  it  suffices  to  show  that  ka 
as  defined  in  (3.3.7)  is  square  integrable  on  Ml2  Now,  with  a  formal  change  in  the 
order  of  integration  and  the  use  of  Parseval’s  identity  the  far  right-hand-side  of 
(3.3.8)  can  be  recast  as 

+00  +00  +  00  +  00 

k2a(x2,rj)  dr]  dx2  = 

—  oo  — oo  — oo  — oo 

+  oo  4*  oo 

=  2~  J  J  \^{k»Ku^v)\ 2dudq.  (3.3.9) 

— oo  — oo 


// 


k](x2,r))dx2  dr] 


The  function  ^r{^s}(-,?7)  :  1R —*  (T  which  appears  in  (3.3.9)  represents,  for  each  tj 
in  Ml,  the  Fourier  transform  of  ka(-,rj).  This  is  supplied  by 

+  oo 

F{ka}(u>,  r])=  J ka(x2,r])e~lu,X2  dx2 

—  OO 

=s(u)1  +  tU}TI-~-—  V(w,i ;)€«*  (3.3.10) 

V 
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where  s  :  2R  — >  (F,  in  turn,  is  the  Fourier  transform  of  s: 


+oo 

s(u)  =  J s(x 2)e~,u,X2  dx 2  Vu;  G  2FZ. 

—  OO 


(3.3.11) 


A  formal  change  in  the  order  of  integration  on  the  far  right-hand-side  of  (3.3.9) 
yields 


+  oo  +oo 


+  oo  +oo 


J  J k](x2,r))dT)dx2  =  ~  J  J  \7{ka}{u,r])?  drjdw, 

—  OO—OO  —00—00 

so  that,  with  the  aid  of  (3.3.10)  and  (3.3.11), 


+  oo  +oo 


+oo  +oo 


J  J k23(x2,T])dT]dX2  =  J  |s(o?)|2(  J 


|1  +  it]w  -  e^wj 2 


27T7;4 

—  OO  — OO 

+oo  -f  oo 


dr])  do 


=  J  |w|3|i(u>)|2  dw  J 


|1  +  i(  -  e'c|2 
27T  C4 


dC  (3.3.12) 


Next,  a  straightforward  application  of  contour  integration  yields  the  identity 


+  °o 

/ 


|l  +  »C-e«|a  ..  1 

2>r<<  '  3’ 


so  that  (3.3.12)  implies  that 


+00  +  00  +oo 

J  J k23(x2,v)dTldx2  =  |  J  M3|s(w)|2<k>.  (3.3.13) 

—  00  —00  — OO 

Note  that,  provided  the  integral  on  the  right-hand-side  of  (3.3.13)  exists,  the 
two  formal  changes  in  the  order  of  integration  performed  above  are  justified  by 
Fubini’s  theorem.5  Now,  by  (3.3.12),  elementary  identities  involving  the  Fourier 

5  See  Halmos  [17]  for  a  statement  and  proof  of  Fubini’s  theorem. 


-34- 


transform  of  first  and  second  derivatives,  and  Parseval’s  identity,  (3.3.13)  gives 

+00+00  +oo 

J  J k23(x2,T))dr)dx2  -  ~  J |s'(x2)||s"(x2)|  dx2.  (3.3.14) 

—  OO  — OO  — OO 

Hence,  provided  s  is  an  element  of  A  it  is  apparent  from  (3.3.8),  (3.3.14)  and  the 
Cauchy-Schwarz  inequality  that  the  kernel  K,  is  square  integrable  on  1R2  and, 
moreover,  that  11^11^2(^2)  can  be  estimated  as  follows: 

-f  OO 

1-^*1  L2(^2)  —  "3”  J \s  (rc2)||-s/,(^2)|  dx 2  <  ^  |-s,||i2(^)  |I5,,|l2(k)  •  (3.3.15) 

—  OO 

Observe  that  while  the  membership  of  s  in  A  is  certainly  sufficient  to  ensure 
that  Ks  is  an  element  of  L2(]R 2)  it  is  not  necessary.  An  application  of  Holder’s 
inequality  to  (3.3.14)  shows,  for  instance,  that  in  order  for  Ks  to  be  square  inte¬ 
grable  on  M 2  it  is  sufficient  to  require  that  s'  and  s"  be  elements  of  Lp( 1R )  and 
Lq(lR),  respectively,  for  some  p  in  (l,oo)  and  conjugate  exponent  q  =  p/(p  —  1). 
The  choice  p  =  q  =  2  clearly  leads  to  the  ultimate  estimate  in  (3.3.14).  It  is, 
moreover,  clear  that,  provided  s  is  in  A,  the  domain  and  range  of  the  operator 
Ms  introduced  in  (3.2.13)  can  both  be  taken  as  L2(JR).  An  immediate  conse¬ 
quence  of  this  observation  is  thai  if  1 p  is  in  L:\JR)  then  so  also  is  M3tP  for  any 
natural  number  n. 

Next,  given  that  Ka  is  an  element  of  L2(1R2)  for  every  function  s  in  A, 
consider  the  issue  of  proving  that  the  forcing  /,  is  similarly  square  integrable  on 
its  domain  of  definition.  Observe,  first,  that  the  singular  behavior  of  Ia  which 
appears  in  the  definition  (3.2.15)  suggests  that  fa  can  be  linearly  decomposed 
into  a  regular  part  and  a  Cauchy  principal  value  part  as  follows: 

+00  +00 

/.(«)»  /(/.to,*)-— J-Tym+ 

J  X2  -  ?  J  X2  ~  £ 

—  00  —00 

+  00  +00 

=  -  /*,(**, /^f 

J  X2  ~  £  J  x 2  -  f 

—  00  —00 


Vz2  €  1R.  (3.3.16) 


-35- 


It  is  now  convenient  to  define  functions  gs  :  1R  — ►  jR  and  hs  :  JR  — »  2R  by 


+  oo 


».(*»)--  I S}()s'(Odi  Vx2£1R, 

J  X2  ~  C 


(3.3.17) 


and 


h„{x2 ) 


+oo 

=  f3'®' 

J  x  2  - 


(o^e 
e 


Vx2  e  1R, 


(3.3.18) 


respectively.  Consider  the  term  of  the  decomposition  involving  the  function  gs. 
From  the  assumed  smoothness  of  s,  the  difference  quotient  which  appears  in  the 
integrand  on  the  right  hand  side  of  (3.3.17)  satisfies 


$(x2)-s(0 


X2  -  £ 


—  M  II Z-00 ( 2R)  ^(x2,0  £  JR  , 


and,  hence,  it  follows  from  the  Cauchy-Schwarz  inequality  that 


+  oo  +00+00 

J di(x2)  dx2  <  iKlloo^)  /</  K»(x2,£)s'(()d02  dx2 


—oo  —  oo 
4-00 -roo 


-f  OO 


/</  K)(xU)iO(  J  lAQpdQdx, 


—oo  —  oo 


5:  I I5,I1>(jr)* 


(3.3.19) 


Therefore,  since  s'  is  contained  in  L2(JR)nL°°(lR ),  (3.3.19)  and  the  bound  (3.3.15) 
on  I*.  I  Li(]R7)  guarantee  that  gs  is  square  integrable  on  2R  and,  furthermore, 
deliver  the  estimate 

1?*Iz,j(r)  -  i’®,lZa(«)  l<s,,lii(jR)  •  (3.3.20) 

Next,  observe  from  (3.3.18)  that  ha  is  a  scalar  multiple  of  the  Hilbert  transform 
of  s'  and  recall  that  the  Hilbert  transform  maps  the  space  Lq(]R)  into  itself  for 
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q  in  (l,cx>).6  Consequently,  the  square  integrability  of  s'  on  Ft  shows  that  ha  is 
also  an  element  of  L2{Ft).  Since  the  space  of  square  integrable  functions  on  Ft 
in  linear  it  is  clear  that  /s,  as  the  sum  of  two  functions  contained  in  L2(Ft),  is 
itself  an  element  of  L2(Ft).  Hence,  under  the  assumption  that  s  is  a  member  of 
the  set  A  the  forcing  fa  must  necessarily  be  square  integrable  on  Ft.  Note  that — 
based  on  the  last  statement  and  earlier  remarks  pertaining  to  the  domain  and 
range  of  the  integral  operator  Ma — the  function  Mafa  is  contained  in  L2(Ft)  for 
every  non-negative  integer  n.  Hence,  if  the  Neumann  series  (3.3.3)  is  uniformly 
convergent  then  <p  must  also  be  square  integrable  on  Ft. 

Up  to  this  point  it  has  been  shown,  as  proposed  above,  that  K,  and  fa  are 
square  integrable  on  their  domains  of  definition  for  every  s  contained  in  A  D  2. 
Finally,  it  is  readily  apparent  from  the  primary  estimate  of  || A.’a  ||  ^ 2 (  jrv j  given  in 
(3.3.15)  that  if  s  is  an  element  of  the  set  X  introduced  in  (3.3.6)  then  inequal¬ 
ity  (3.3.5)  must  hold.  Hence,  (3.3.5)  is  satisfied  for  every  element  s  of  2.  To 
recapitulate,  observe  that  if  s  is  in  1  then  there  exists  a  unique  solution  to  the 
corresponding  integral  equation  (3.2.14)  given  by  the  appropriate  Neumann  series 
(3.3.3). 

There  may  exist  solutions  to  (3.2.15)  which  are  not  obtainable  via  the  Neu¬ 
mann  series  construction.  Since,  however,  the  solution  to  the  integral  equation 
obtained  via  this  construction  is  unique  for  each  s  in  2  it  is  apparent  from  the 
above  discussion  that,  should  there  exist  any  solutions  to  (3.2.15)  which  can  be 
aquired  by  alternate  means,  these  must  correspond  to  curves  C,  described  by 
functions  s  which  do  not  belong  to  J  (and  may  not  even  belong  to  A).  It  is 
interesting  to  speculate  on  whether  some  of  these  solutions  might  correspond  to 
states  wherein  the  phase  boundaries  manifest  large  slopes  and/or  curvatures  akin 
to  those  exhibited  by  the  fingers  found  in  studies  of  porous  media,  solidification, 
and  crystal  growth. 

Prior  to  concluding  this  section  a  few  comments  regarding  the  uniqueness  of 


6 


This  fact  is  established  in  RlESZ  [22], 
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the  solution  to  the  reduced  boundary  value  problem  are  in  order.  Let  s  be  an 
element  of  I.  Then,  by  the  foregoing  results,  the  related  reduced  boundary  value 
problem  has  a  solution  given  by  the  appropriate  Neumann  series  (3.3.3).  It  is 
known  that  the  solution  to  the  integral  equation  which  issues  from  the  reduced 
boundary  value  problem  is  unique.  The  uniqueness  of  the  solution  to  the  reduced 
boundary  value  problem  is,  however,  still  in  question.  It  will  now  be  shown 
that  the  solution  of  the  reduced  boundary  value  problem  is  unique — just  as  with 
the  globally  elliptic  pairwise  homogeneous  equilibrium  states — up  to  an  arbitrary 
additive  constant.  To  see  this  suppose  that  vi  :  ]R2\Ca  — »  JR  and  v2  :  Si2\Ca  — +  2R 
are  both  solutions  to  the  reduced  boundary  value  problem  corresponding  to  s  in 
X;  define  w  :  1R2  \  Ca  — >  1R  by  their  difference  (uj  —  w2 )  on  1R2  \  Ca.  Then,  from 
(3.1.17),  w  clearly  satisfies  the  following  boundary  value  problem: 

w,aa  -  0  on  JR2\C„ 

dw  ,  x  dw .  .  .  v  m 

/ii^-(5(-)+,-)  =  ^2 -X- («(•)->•)  on  ]R, 

0n  °n  (3.3.21) 

u,(a(-)+>’)  =  »(«(•)-,•)  on  1R, 

w,a  (xi,  -)e<*  =  o(l)  as  x\  — »  ±oo  on  ]R. 

From  (3.2.16)  it  is  readily  apparent  that  a  solution  to  (3.3.21)  is  provided,  modulo 
an  arbitrary  additive  constant,  by 

+  oo 

w(x\,X2)  =  -l  V(ii,x2)  G  1R2.  (3.3.22) 

—  OO 

where  G®  :  (JR2  \  Ca )  x  JR  — >  JR  is  given  by  (3.2.3)i,  and  %l>  :  JR  — >  JR  satisfies 

0  +  A  Maip  =  0  on  JR.  (3.3.23) 

It  is  clear,  based  on  the  assumption  that  s  is  in  J,  that  A  cannot  be  an  eigenvalue 
of  the  operator  Ma.  Hence,  (3.3.23)  has  only  the  zero  solution.  Now,  since 
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the  representation  (3.3.22)  for  w  is  modulo  an  arbitrary  additive  constant,  the 
functions  iq  and  v-i  can  differ  at  most,  as  stated  above,  by  a  constant. 

During  the  course  of  this  section  it  has  been  shown  that,  for  each  function  s 
contained  in  X  there  exists,  up  to  an  arbitrary  additive  constant,  a  unique  solution 
to  the  reduced  boundary  value  problem  (3.1.17).  This  solution  corresponds  to 
a  deformation  involving  a  non-planar  interface  Qs.  Inasmuch  as  the  reduced 
phase  segregation  requirement  (3.1.18)  has  not  yet  been  applied  it  is  still  unclear 
whether  any  of  the  aforementioned  solutions  give  rise  to  globally  elliptic  two-phase 
equilibrium  states.  The  next  section  will,  therefore,  focus  on  characterizing  a 
subset  of  X  for  which  there  exist  solutions  to  the  reduced  boundary  value  problem 
augmented  by  the  (reduced)  constraints  of  phase  segregation.  If  a  function  s 
belongs  to  this  subset  of  X  the  interface  Qs  will  qualify  as  a  phase  boundary. 

3.4.  Implementation  and  satisfaction  of  the  reduced  phase  seg¬ 
regation  requirement.  Let  s  be  an  element  of  X  and  suppose  that  and 
(up  to  an  additive  constant)  v  are  the  corresponding  solutions  to  the  integral 
equation  (3.2.15)  and  the  reduced  boundary  value  problem  (3.1.17).  If  v  is  to 
provide — through  (3.1.16) — a  solution  u  to  the  constrained  boundary  value  prob¬ 
lem  its  gradient  must  comply  with  the  reduced  strain  constraints  (3.1.18).  Let 
k  :  R 2  \  Ca  — ♦  R+  denote  the  reduced  shear  strain  field  given  by 

k  =  v,qv,q  on  R2  \C».  (3.4.1) 

Certainly  |u,i|  must  be  less  than  or  equal  to  k  on  R2  \  Cs;  hence,  if  the 
reduced  shear  strain  field  complies  with 

<  min{7  -  yr,7i  -  7}  on  R2  \CS,  (3.4.2) 

* 

then  both  of  the  inequalities  which  comprise  the  reduced  phase  segregation  re¬ 
quirement  (3.1.18)  will  be  satisfied.  Notice  that  the  foregoing  condition  is  suffi¬ 
cient  but  not  necessary  to  ensure  the  segregation  of  phases.  It  may,  consequently, 
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► 


lead  to  overly  conservative  restrictions.  Despite  the  strong  restrictions  which  may 
be  imposed  by  enforcing  (3.4.2)  in  lieu  of  (3.1.18),  it  will  be  demonstrated  that 
there  exists  a  non-empty  subset  of  J  each  element  of  which  gives  rise  to  a  soluble 
reduced  constrained  boundary  value  problem  with  a  reduced  shear  strain  field  k 
that  allows  their  satisfaction. 

The  following  simple  calculation  shows  that  k  is  subharmonic  on  JR2  \  C, : 

^iaa  =  (^)/9  ))oa  —  2(u,ayj  W,/3  ))Q, 


—  2 (v,o/J  Via/}  +V,QQp  V,0  ) 

—  2u,q^  ^  0  on  JR  \  Cj. 


The  subharmonicity  of  k  on  JR2  \  Cs  implies,  given  the  decay  properties  of  the 

O  .  O 

gradient  of  v  embodied  by  (3.1.17)4,  that  its  maximum  values  on  and  T>1 
must  occur  in  the  limits  approaching  the  curve  Cs  from  the  high  and  low  strain 
sides,  respectively.  Hence,  in  determining  whether  the  reduced  shear  strain  field 
k  satisfies  (3.4.2)  it  is  sufficient  to  analyze  its  limiting  behavior  on  either  side  of 
the  curve  C„.  A  convenient  approach  to  this  is  afforded  by  examining  the  limits 
of  the  normal  and  tangential  derivatives  of  v  on  either  side  of  Cs.  FVom  (3.2.6), 
(3.2.11)  and  (3.2.12)  it  is  evident  that  the  limiting  values  of  the  normal  derivative 
of  the  the  reduced  out-of- plane  displacement  field  are  given  by 


dv 

dn 


(s(x2)-,x2) 


on  the  high  strain  side  of  Cs,  and 


7/y(*2) 

y/l  +  S'(X  2)2 


Vx2  €  JR 


(3.4.3) 


dv 

dn 


(s(x2)+,x2) 


lMx2) 

y/l  +  s'(x2)2 


Vx2  €  JR 


(3.4.4) 


on  the  low  strain  side  of  Ct. 

Let  l :  JR.  Af  designate  the  unit  tangent  vector  to  Qs  defined  by 


l(*a)  = 


s,(j2)ei  +  e2 
y/\  +  s'(x2)2 


Vx2  €  JR. 


(3.4.5) 
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Then  a  calculation  very  similar  to  that  used  in  obtaining  (3.2.4)  yields 


dv 

Tl 


(s(x2)±,X2)  = 


7/  -  7r 

2  Ks/\  +  s'(x2)2 


+oo 

— oo 


+ 


7t  ~  7r 

27rv/l  +  s'(x2)2 


+oo 


±  — .2L  =  7r^z::s/(x2)  Vx2  €  2R  (3.4.6) 

2y/l  +  s'(x2)2  ^ 


for  the  limits  of  the  tangential  derivative  of  u  on  either  side  of  the  curve  Cs.  Here, 
as  in  (3.2.4),  /,(x2,  •)  and  ifj(x2,  •)  are  given,  for  each  x2  in  IR,  by  (3.2.5). 

Turn  now  to  the  estimation  of  (3.4.3),  (3.4.4)  and  (3.4. 6)±.  Consider  the 
limits  of  the  normal  derivative  first.  The  following  pair  of  inequalities  follow 
immediately  from  (3.4.3)  and  (3.4.4): 


dv 

dn 


(s(x2)-,x2) 


<7/|v?(x2)|  <  7/  Hv?lljL~(iR) 


Vx2  G  2R, 


dv 

dn 


(s(x2)+,x2) 


<  7rM®2>|  <  7r|v?lL»(JR) 


Vx2  e  m. 


(3.4.7) 


Hence,  in  order  to  bound  the  limits  of  the  normal  derivative  of  the  reduced  dis¬ 
placement  field  on  either  side  of  Ca  it  is  only  necessary  to  estimate  the  L°°  norm 
of  tp  over  1R.  In  Appendix  A  it  is  shown  that  if  0  :  2R  — ►  iR  is  an  element  of  the 
set  V  defined  in  (3.1.7)  then  ||0l£oo(j*)  exists  and  can  be  bounded  as  follows: 


101l«>(K)  ^  2(-|)4  ||V’|£s(jr)  10  Il*(ir)  •  (3.4.8) 


Recall  from  Section  3.3  that  if  s  is  contained  in  T  then  <p  is  square  integrable 
on  IR.  Therefore,  if  <p'  exists  and  is  square  integrable  on  2R  the  inequality  displayed 
in  (3.4.8)  can  be  used — with  0  replaced  by  <p — to  obtain  an  estimate  for  the  L°° 
norm  of  ip  over  M.  Suppose,  from  now  on,  that  s  is  three  times  continuously 
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differentiable  element  of  T  with  a  square  integrable  third  derivative  on  R.  It 
can  be  readily  shown  that  this  is  sufficient  to  guarantee  that  tp'  exists  and  is  an 
element  of  L2(R).  In  Appendix  B  it  is  demonstrated  that  the  L2  norms  of  <p  and 
<p'  over  R  can  be  estimated,  respectively,  by 


lull , 


Cl(l 


(3.4.9) 


and 


R)  -  c2  [  Is"Il*(1R) 

+  Is  IIl2(JR) 

+  II6  lia(jR)  \s  ||£2(iR)  Mll,*(]R) 

+  IIs  IIl2(2R)  l's//|lz,*(2R)  Is,,,|Il*(JR)  IMIl*(.R)]-  (3.4.10) 


The  constants  ci  and  c-j.  which  appear  in  (3.4.9)  and  (3.4.10)  are  positive  real 
numbers  entirely  independent  of  s.  Note  that  the  denominator  in  (3.4.9)  is  strictly 
positive  since  s  is  an  element  of  J.  With  the  aid  of  (3.4.9)  the  estimate  (3.4.10) 
for  can  be  expressed  completely  in  terms  of  |-s# )Iz,2(ir)»  I5”  I  £,*(«)  and 

II5,/,|Il*(jr)'  Hence,  (3.4.9),  (3.4.10),  (3.4.8)  and  (3.4.7)  give  estimates  for  the 
moduli  of  the  limiting  values  of  the  normal  component  of  the  gradient  of  v  on 
either  side  of  C9  in  terms  of  the  L2  norms  of  the  first  three  derivatives  of  s  over 
R. 

To  provide  an  estimate  for  it  it  ramains  to  obtain  bounds  on  the  limiting 
values  of  the  tangential  derivative  of  v.  Toward  this  objective,  introduce  a  function 
A  :  R  — ►  R+  by 


+  oo  +oo 

A(x2)  =  J K.(x2,()s'( ()<%  +  j  /,(*2,£MO^  +  ™'(*2)  Vx2  G  R.  (3.4.11) 

— oo  — oo 
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Then  it  is  clear  from  (3.4.6)  and  (3.4.11)  that 

<  Vi2  6  M.  (3.4.12) 

Under  the  current  assumption  that  s'"  exists  and  is  square  integrable  on  JR 
it  can  be  shown  that  both  A  and  A'  are  elements  of  L2(]R).  Inequality  (3.4.8) 
can,  therefore,  be  applied  with  ip  replaced  by  A ;  this  leads,  through  (3.4.12),  to 
a  bound  on  the  limiting  values  of  the  tangential  component  of  the  gradient  of 
v  on  either  side  of  Ca.  Given  the  bounds  (3.4.9)  and  (3.4.10)  for  and 

II^1l2(.R)  ^  is  straightforward  to  derive  estimates  for  the  L2  norms  of  A  and  A! 
over  JR  in  the  form 

Mll,2(2R)  —  c3  [  I5  Il2(2R)  +  II5  Il2(2R)  I5  H2(K)  +  Iv^l L2(2R) 

+  I5,|Il2(K)  l^lL2(iR)]>  (3.4.13) 


dv 

x2)±,x2) 


and 


1^  i jL2 ( 2R)  —  C4  [l5  li2(iR)  +  II5  Il2(J{)  IIs  Il2(«)  II5  l£2(JR) 

+  II5  ll£2(H)  I-  Ilj (2R)  ll5”  I£2(k) 

+  ll5  li2(iR)  II5  II Z/2 (2R)  l5”,|l£2(JR)  IMIl2(JR) 

+  Is  II£2(k)  Is  IIl2(2R)  is,”l£2(ji)  I¥,IIl2(k) 

+  IMIl2(jr)]-  (3.4.14) 

respectively.  Here  C3  and  c4  are  positive  real  numbers  which,  like  c\  and  C2,  are 
independent  of  s.  Note  that  the  L 2  norm  of  9?  over  JR  appears  in  the  estimates 
for  both  \A\L2(R}  and  Hence,  with  the  aid  of  (3.4.9)  these  can  be 

rewritten  solely  in  terms  of  the  L 2  norms  of  s' ,  s"  and  s'"  over  JR.  As  with 
the  normal  component  of  the  gradient  of  v ,  a  bound  on  the  limiting  values  of 
the  tangential  component  of  the  gradient  of  v  on  either  side  of  Cs  is  given,  from 
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(3.4.13),  (3.4.14),  (3.4.11),  (3.4.12),  (3.4.9)  and  (3.4.8)  in  terms  of  the  L2  norms 
of  the  first  three  derivatives  of  s  over  JR. 

From  the  proceeding  discussion  it  is  clear  that  k  can  be  made  arbitrarily  small 
by  reducing  the  size  of  |s,||i2(JR),  |s,#Il*(.h)  and  More  specifically, 

(3.4.9),  (3.4.10),  (3.4.13)  and  (3.4.14)  can  be  substituted  appropriately  into  (3.4.8) 
to  compute,  using  (3.4.7)  and  (3.4.12),  an  upper  bound  for  k  on  £>;  U  T>1  in  the 
form 

K  <  ■F2(|s'||i2(jR^ ,  |s"|L2(^) ,  =‘.r2.  (3.4.15) 

Define  a  set  of  functions  J  by 

J  =  {-s  €  -4 1  ra  <  min {7  -  7,., 7/  -  7}}  n  X,  (3.4.16) 

* 

where  X  is  given  by 


X  =  {x :  R-*  JR|V>  €  C3(JR),i/>(n)  €  L2(JR),n  =  0,1, 2, 3).  (3.4.17) 

Then,  from  (3.4.2),  (3.4.15)  and  (3.4.16),  provided  s  is  an  element  of  the  set  77 
defined  by 

i7  =  lnj,  (3.4.18) 

there  exists  a  solution — unique  up  to  an  additive  constant — to  the  associated  con¬ 
strained  boundary  value  problem  (3.1.17)-(3.1.18)  with  phase  boundary  Qs.  This 
solution  defines  a  globally  elliptic  inhomogeneous  two-phase  equilibrium  state  and, 
therefore,  establishes  the  existence  result  sought  after  here.  Note — from  the  defi¬ 
nition  of  II — that  the  approach  delineated  above  provides  a  means  by  which  an 
uncountably  infinite  number  of  such  states  can  be  constructed.  It  is  significant 
that  the  loading  conditions  related  at  the  outset  of  Section  3.1  give  rise  to  not 
only  a  globally  elliptic  pairwise  homogeneous  equilibrium  state  but  also  an  un¬ 
countably  infinite  number  of  globally  elliptic  inhomogeneous  two-phase  equilibria. 
This  result  clearly  reflects  the  underlying  non-linearity  of  the  problem. 
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As  remarked  earlier,  there  may  exist  globally  elliptic  inhomogeneous  two- 
phase  equilibrium  states  which  cannot  be  constructed  via  the  approach  taken 
above  and,  thus,  do  not  correspond  to  phase  boundaries  in  the  set  i7.  Under 
relaxed  smoothness  assumptions  on  s  there  may,  however,  exist  still  other  globally 
elliptic  inhomogeneous  two-phase  equilibrium  states  which  can  be  constructed 
via  Neumann  series.  In  particular,  under  such  relaxed  circumstances,  it  may 
be  possible  to  demonstrate  the  existence  of  states  wherein  the  associated  phase 
boundaries  exhibit  geometrical  irregularities  such  as  corners  or  cusps  (recall  that 
the  existence  of  equilibria  involving  cusped  phase  boundaries  has  been  established 
by  Rosakis  [23]  in  his  work  involving  a  special  anisotropic  material). 

3.5.  An  example.  Given  the  results  of  Sections  3.3  and  3.4  it  is  illuminating 
to  consider  a  particular  class  of  functions  in  the  set  A  defined  by  (3.1.6)  and 
determine  a  subset  of  this  class  of  functions  which  are  also  contained  in  the  set 
II  defined  by  (3.4.18).  Toward  this  end,  let  s  is  given  by 

s(x 2)  =  1  Vxz  6  R'  (3.5.1) 

where  h  and  i  are  both  positive  constants.  A  representative  graph  of  s  is  displayed 
in  Figure  3.  Note  that  s  is  clearly  an  infinitely  differentiable  element  of  the  set  A 
for  all  values  of  the  parameters  h  and  l.  Let  the  ratio  of  h  to  l  be  denoted  by  e. 
The  kernel  Ka  associated  with  s  must,  as  a  consequence  of  the  results  of  Section 
3.3,  be  square  integrable  on  1R2.  In  fact,  from  (3.3.15)  one  finds,  after  a  bit  of 
calculation,  that  the  L 2  norm  of  Ka  over  M2  can  be  bounded  as  follows: 

=  2e‘  (3.5.2) 

In  the  latter,  ka  is  as  defined  in  (3.3.7).  Hence,  it  is  clear  from  (3.2.10) 
and  (3.3.6)  that  if — for  a  given  choice  of  the  moduli  Hi  and  /x2  which  define  the 
elliptic  phases  of  the  three-phase  material  with  shear  stress  response  function 


Tp — the  parameter  e  satisfies 


c  <  2  --  /i2 ,  (3.5.3) 

^1-/^2 

then  the  function  s  introduced  above  will  be  an  element  of  X  fl  C°°{IR).  Assume, 
from  now  on,  that  s  as  defined  in  (3.5.1)  is  such  that  (3.5.3)  holds.  Then  the 
Neumann  series  (3.3.3)  converges  uniformly  on  U  to  a  solution  of  the  intergal 
equation.  With  the  aid  of  the  decomposition  (3.3.16)  of  the  forcing  /,,  the  solution 
of  the  integral  equation  can  be  expressed  as 

OO  OO 

^(z2)  =  A^(-A)"(MJn/lj)(2;2)  +  A^(-A)n(MJny3)(a:2)  Vx2  G  JR,  (3.5.4) 

n=0  n=0 


where  M.a  is  as  defined  in  (3.2.13).  Given  rp  and,  hence,  the  moduli  /ii  and  //2,  it 
can  be  readily  shown  that,  for  every  e  which  satisfies  (3.5.3),  the  following  order 
relations  hold  for  each  non-negative  integer  n: 


(Mn3ha)(x2)  =  0(e2n+1)  Vx2  G  JR, 
(MMix 2)  =  0(£2n+2)  Vx2  €  JR, 


(3.5.5) 


Therefore,  facilitated  by  (3.5.4)  and  (3.5.5),  <p  can  be  represented  in  the  form 


•f  OO 

Y>(x2)  =  2A hi2  /  - 


(,di 


{p  +  em_X2) 


-f-  0(e2)  Vx2gJR. 


(3.5.6) 


An  application  of  contour  integration  yields 


+oo 

/ 


v2 + ent  -  *2) 


n  l~(^)2 

2^3(l  +  (^)2)2 


Vx2  G  JR, 


so  that  (3.5.6)  becomes,  with  the  aid  of  (3.2.11), 


V?(x2)  = 


//l  -M2  l-(¥)2 

+^2  (H-(^)2)2 


e  +  0(e2) 


Vx2  €  JR. 


(3.5.7) 
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Observe  that  (3.5.7)  and  (3.5.1)  can  be  used  in  (3.2.16)  and  (3.1.16)  to  con¬ 
struct,  for  each  appropriate  pair  (h,£),  an  approximate  solution  to  the  reduced 
boundary  value  problem  (3.1.17).  Now,  substitution  of  (3.5.1)  and  (3.5.7)  into 
(3.4.3)  and  (3.4.4)  delivers  the  following  formulae  for  the  limiting  values  of  the 
normal  derivative  of  v  on  either  side  of  Ca: 

£  «»»>-■*»>  -  (i+(ff)* £  +  0(e2)  VX2£K' 


(3.5.8) 

#.,+*  (i+(W£+0(£2)  Vl2eJR- 

Similarly,  substitution  of  (3.5.1)  and  (3.5.7)  into  (3.4.6)  gives  rise  to  the  following 
expressions  for  the  limiting  values  of  the  tangential  derivative  of  v  on  either  side 
of  C3: 

2/xi(7j  —  7r)  2n  w  ^ 

+  (.£1)2)2  e  +  °(e  )  Vx2  G  JR, 

(3.5.9) 

e  +  0(e2)  Vx2  6  JR. 


jf(s(x2)-,x2)  = 


^(s(x2)+,  x2)  = 


Hi  +  Hi  (1  +  (^f)2)2 
2/X2(7l  -  7r) 


/ 


/il+M2  (l  +  (-^)2)2 

The  expansions  in  (3.5.8)  and  (3.5.9)  show  the  dependence  on  e  of  the  limiting 
values  of  the  normal  and  tangential  derivatives  of  v  on  either  side  of  Cs.  They 
readily  imply  that,  for  the  function  s  indicated  in  (3.5.1),  the  associated  reduced 
shear  strain  field  k  introduced  via  (3.4.1)  satisfies 


«  =  0(e2)  on  JR\CS 


(3.5.10) 


An  immediate  consequence  of  (3.5.10)  is  that  if  e  is  made  sufficiently  small  the 
reduced  phase  segregation  requirement  (3.1.18)  will  then  be  satisfied  and  the 
relevant  function  s  will  be  contained  in  J7. 

Note,  alternatively,  that  the  L 2  norms,  over  JR,  of  the  first  three  derivatives 
of  the  function  s  defined  in  (3.5.1)  can  be  computed  directly  to  give 


I5"Nl*(«)  - 


Is'"  II 


L*(R)  ~ 
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If  the  foregoing  axe  substituted  in  the  estimates  (3.4.9),  (3.4.10),  (3.4.13),  and 
(3.4.14)  then  it  is  straightforward  to  show,  with  appropriate  use  of  (3.4.7),  that  the 
quantity  T,  defined  in  (3.4.15)  is  of  order  e — which  corroborates  the  asymptotic 
results  obtained  above.  Hence,  if  the  parameters  h  and  I  which  appear  in  the 
definition  of  s  are  chosen  so  that  e  is  sufficiently  small,  f,  will  satisfy 


r3  <  min{7  -  7r,  7/  -  7  }» 


(3.5.11) 


and  s  will  be  an  element  of  17. 

In  either  case  a  class  of  phase  boundaries  Qs  for  which  e  =  h/i  is  sufficiently 
small  emerges  from  the  class  of  functions  s  given  by  (3.5.1).  The  reduced  out- 
of-plane  displacement  field  corresponding  to  each  such  s,  is  by  using  (3.5.1)  and 
(3.5.4)-(3.5.7)  in  (3.2.16),  given  approximately  by 


•u(xi ,  X2)  ~ 


7/  -  7r  hx2 


+00 

/ 


dt 


2*  i  J  a +e2x(t)2+(t-r) 
i-e2 


—00 
+00 


V  -It 
27 r 


/ 


(i+e2)2 


In 


V(*i,x2)e2R2\C,.  (3.5.12) 


An  approximation  to  the  corresponding  primitive  out-of-plane  displacement  field 
u  is  then  calculated  easily  by  substituting  (3.5.12)  appropriately  into  (3.1.16). 
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4.  STUDY  OF  PHASE  BOUNDARY  KINETICS  AND  STABILITY 

This  chapter  relies  on  the  concept  of  a  quasistatic  motion  introduced  in  Sec¬ 
tion  2.3.  Recall  that,  in  addition  to  the  shear  stress  response  function  rp,  the 
constitutive  characterization  of  the  material  at  hand  includes  a  kinetic  response 
function  V  which,  in  the  setting  of  a  quasistatic  motion,  dictates  the  dependence 
of  the  normal  velocity  of  a  particle  located  at  a  point  on  a  phase  boundary  on  the 
driving  traction  acting  at  that  point.  Given  the  distribution  of  driving  traction  on 
a  particular  phase  boundary  it  is  therefore — through  the  kinetic  relation — possible 
to  discuss  the  kinetics  and  stability  of  that  phase  boundary  in  slow  motions.  For 
illustrative  purposes  this  is  done  below  in  the  context  of  the  specific  class  of 
phase  boundaries  studied  in  Section  3.5.  A  similar  analysis  could,  in  principle,  be 
performed  for  any  function  s  in  77. 

In  Section  4.1  the  driving  traction  /  which  acts  on  such  a  phase  boundary  is 
derived.  It  is  demonstrated  that  /  is  composed  of  the  sum  of  an  ambient  term  fo 
which  corresponds  to  the  constant  driving  traction  which  would  act  on  a  planar 
phase  boundary  corresponding  to  a  suitable  globally  elliptic  pairwise  homogeneous 
equilibrium  state  and  higher  order  terms  which  represent  the  increment  to  the 
driving  traction  resulting  from  the  non-planarity  of  the  surface  Qs. 

In  Section  4.2  the  ambient  term  fo  and  the  most  significant  non-constant 
term  in  the  driving  traction  /  are  used  in  conjunction  with  V  and  the  kinetic 
relation  to  address  phase  boundary  kinetics  and  stability. 

4.1.  The  driving  traction  acting  on  an  arbitrary  element  of  a  spe¬ 
cific  class  of  phase  boundaries.  Let  s  be  given  by  (3.5.1)  with  e  =  h/t  chosen 
so  that  (3.5.11)  is  fulfilled.  As  discussed  in  Section  3.5,  Qa  is  then  a  phase  bound¬ 
ary.  Consider  the  computation  of  the  driving  traction  acting  on  Qs.  The  simple 
manner  in  which  Qs  can  be  parameterized  implies  that — in  the  present  context 
of  antiplane  shear — the  expression  for  the  driving  traction  provided  in  (2.4.6)  can 
be  written  as  a  function  of  one  variable.  Furthermore,  it  can  be  shown  with¬ 
out  difficulty  that,  for  the  special  three-phase  material  with  shear  stress  response 
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function  rp,  the  driving  traction  /  is  given  by 


/  = 


Ml  ~  M2 


(u,a  ($(•)+, -77)  on  R. 


(4.1.1) 


The  definition  of  v  supplied  in  (3.1.16)  readily  furnishes  the  following  expres¬ 
sions  for  u,Q  (s(-)±,  -)ea  on 


u,a  (s(-)-,  -)e0  =  7/ei  +  v,a  ($(•)-,  -)ec 


=  w +  £«•)-,  •)n  +  fw-)-.-)|  on  m , 


U,c  (s(-)  +  ,  -)ea  =  7rei  +  V,a  (s(-)  +  ,  0ea 


(4.1.2) 


=  7r«i  +  |^(s(-)+?-)n  + ^(5(‘)+r)l  on  R. 


It  is  possible  to  show,  from  (3.5.1),  (3.1.7)  and  (3.4.5),  that  the  unit  normal  and 
tangent  vectors  to  Q„  satisfy  the  following  order  relations: 


n  •  ei  =  1 0(c)  on  R,  1  •  ei  =  0(e)  on  R.  (4.1.3) 


Hence,  (3.1.16),  (4.1.2)  and  (4.1.3)  yield 


dv  dv 

U,a  (s(-)+,  ■)“»«*  («(•)-.  •)  =  7/7r  +  7l^(<(-)+»  •)  +  ') 


+  0(e2)  on  R. 


(4.1.4) 


Now,  if  (3.5.8)  is  inserted  appropriately  in  (4.1.4)  and  the  result  is  substituted 
into  (4.1.1)  the  driving  traction  along  the  phase  boundary  Qs  can  be  expanded 
in  powers  of  c  as  follows 


/(*»)  =  h  +  Vx,€J R, 


(l+(¥)2)2 

with  the  constants  /0  and  v  given  by 


(4.1.5) 


fo  = 


Hi  ~  H2 


(7i7r  -  77).  v  = 


(mi  -  M2)2 
(Ml  +  M2) 


7/7r. 


(4.1.6) 
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/o  is  the  ambient  or  base  driving  traction  which  would  hold  on  a  planar  phase 
boundary  associated  with  a  globally  elliptic  pairwise  homogeneous  equilibrium 
state  with  displacement  field  (3.1.1).  Observe  that  since 


8  -  (l  +  (f  )2)2 


<  1  Vx2  €  R, 


the  coefficient  of  the  0(e)  term  in  the  expansion  of  /  provides  a  bounded  correction 
to  the  base  term  fo ■  This  term  is,  for  small  e,  the  most  significant  contribution 
to  /  which  results  from  the  deviation  in  the  geometry  of  Q„  from  planar.  Note 
that  fo  can  take  on  any  real  value  whereas  v  must  be  positive.  Furthermore, 
because  the  difference  (/ii  —  /z2)  is  squared  in  (4.1.6)2,  the  positivity  of  v  holds 
even  if  the  affiliations  of  the  moduli  and  /x2  are  reversed  so  as  to  be  associated 
with  the  high  and  low  strain  phases  of  the  material  at  hand.  From  (4.1.5)  it  is 
apparent  that  the  0(e)  contribution  to  the  distribution  of  driving  traction  on  the 
phase  boundary  Q„  under  consideration  is  a  complicated  function  of  position.. 
See  Figure  4  for  the  graph  of  the  function  corresponding  to  the  0(e)  term  in  the 
expansion  of  /.  The  results  of  ABEYARATNE  [1],  imply  that  the  equilibria  at  hand 
must,  in  general,  constitute  metastable  states. 

4.2.  Kinetics  and  stability  of  an  arbitrary  element  of  a  particular 
class  of  phase  boundaries.  Given  the  expansion  (4.1.5)  consider  the  issue  of 
analyzing  the  kinetics  and  stability  of  an  arbitrary  element  of  the  class  of  phase 
boundaries  at  hand.  Suppose  that  V  is  twice  continuously  differentiable  on  R. 
Then,  as  noted  in  the  remarks  following  (2.3.9),  the  admissibility  of  V  requires 
that  V’(O)  =  0  and  F'(0)  >  0.  Assume,  for  the  purposes  of  this  discussion,  that  fo 
is  not  a  critical  point  of  V ;  note  that  this  requires,  in  particular,  that  if  fo  =  0  then 
F'(0)  >  0.  See  Figure  5  for  examples  of  graphs  of  monotone  and  non-monotone 
kinetic  response  functions. 

Let  Qcg  and  be  those  subsets  of  the  phase  boundary  Qs  defined  as  follows: 


a;  =  {x€e.iu q{  =  q.\qi 
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Observe,  from  (4.1.5),  that  Qcs  and  Q{  correspond  to  the  portions  of  Qs  upon 
which  the  0(e)  correction  to  /o  is  negative  and  positive,  respectively.  Note, 
also,  that  is  the  subset  of  Q,  whose  geometry  deviates  most  significantly 
from  planar — that  is,  rougly  speaking,  the  major  portion  of  the  bump  which  is 
associated  with  the  graph  of  the  function  s  given  by  (3.5.1)  corresponds  to  the 
image  of  (—  £,£)  under  s.  See  Figure  4. 

From  the  assumed  smoothness  of  V,  (4.1.5)  and  Taylor’s  theorem  the  normal 
velocity  at  a  point  on  the  phase  boundary  is  given  by 

V.W  =  ?(/(»»))  =  V(M  -  -V'(/o)^~((g)2  e  +  0(S)  Vx6S,.  (4.2.1) 

In  determining  the  kinetic  tendencies  of  Q,  it  is  now  convenient  to  consider  two 
cases.  These  are  /o  =  0  and  fo  ^  0.  Note  that  the  base  globally  elliptic  pairwise 
homogeneous  equilibrium  state  is  mechanically  stable  only  in  the  first  of  these 
two  cases. 

Consider  the  case  fo  =  0.  Then,  since  V"(0)  =  0,  (4.2.1)  imples  that 

V-„(x)  =  t  +  0(e2)  Vx  e  Q„  (4.2.2) 

Since  V'(0)  and  v  are  positive,  it  is  apparent  from  (4.2.2)  and  (4.1.3)i  that,  to 
most  significant  order  in  e,  all  points  on  Ql  tend  to  move  in  the  — ei  direction 
while  all  points  on  Q{  tend  to  move  in  the  ei  direction.  That  is,  if  f0  =  0  then 
the  phase  boundary  displays  a  proclivity  to  become  planar. 

Now  consider  the  case  where  /o  ^  0.  Suppose,  first,  that  /o  is  positive.  As 
such  the  dominant  contribution  to  the  normal  velocity  is,  at  all  points  on  Qa ,  in 
the  ej  direction.  Recall  that  /o  is  assumed  not  to  be  a  critical  point  of  V ;  hence, 
since  /o  ^  0,  V'(fo)  can  be  either  positive  or  negative.  If  1^(0)  >  0  then,  since 
v  >  0,  the  normal  velocity  of  points  on  Qct  and  Q{  will  decrease  and  increase, 
respectively,  on  top  of  the  ambient  value  V(fo)-  This,  as  in  the  case  where 
fo  =  0,  indicates  a  tendency  for  the  phase  boundary  to  straighten  out.  If,  however, 


-52- 


V'(fo)  <  0  then,  since  v  >  0,  the  exact  opposite  occurs — the  normal  velocity  of 
points  on  Qcs  and  Q{  will  add  positive  and  negative  increments,  respectively,  to  the 
ambient  value  V(fo).  The  protruding  part  of  the  phase  boundary,  if  V'(fo)  <  0, 
portrays  a  tendency  to  grow  while  the  flat  part  lags  behind.  The  subcase  where 
/o  is  negative  is  yields  a  completely  analogous  result.  That  is,  when  /o  <  0 
the  phase  boundary  shows  a  propensity  to  become  planar  or  develop  a  larger 
protrusion  depending  upon  whether  V'(fo)  is  positive  or  negative,  respectively. 

The  foregoing  discussion  shows  that  the  kinetics  of  a  phase  boundary  Qt 
in  the  class  at  hand  are,  to  first  order  in  e,  stable  or  unstable  depending  upon 
whether  the  kinetic  response  function  is  locally  increasing  or  decreasing  at  the 
ambient  driving  traction  /o.  If  the  constitutive  description  of  a  three-phase  ma¬ 
terial  with  shear  stress  response  function  rp  also  includes  a  monotone  increasing 
kinetic  response  function  such  as  that  depicted  in  Figure  5a  it  is  clear  that  phase 
boundaries  of  the  class  under  consideration  will  always  be  stable.  If,  on  the  other 
hand,  the  constitutive  description  incrudes  a  non-monotone  kinetic  response  func¬ 
tion  like  that  depicted  in  Figure  5b  it  is  always  possible  to  choose  7/  and  7r  so 
that  the  phase  boundary  is  unstable.  These  results  suggest  that  it  may  be  reason¬ 
able  to  classify  those  three-phase  materials  with  shear  stress  response  function  rp 
as  kinetically  stable  and  kinetically  unstable  depending  on  whether  the  kinetic  re¬ 
sponse  function  V  is  a  monotone  or  non- monotone  function  of  its  argument.  Such 
a  classification  is  consistent  with  that  found  by  Fried  [13]  in  a  linear  stability 
analysis  of  planar  phase  boundaries  in  arbitrary  three-phase  materials  subjected 
to  a  class  of  perturbations  which  encompasses  the  set  of  phase  boundaries  il 
determined  in  Chapter  3. 
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APPENDICES 


Appendix  A.  In  this  appendix  inequality  (3.4.8)  is  established  for  all  func¬ 
tions  t/>  contained  in  the  set  V  defined  in  (3.1.7).  Let  x  :  JR  — > ►  JR  be  an  element  of 
V  with  compact  support  about  the  origin;  suppose,  further,  that  x(0)  =  1.  Then, 
if  ip  is  contained  in  V,  one  has  the  following  inequality:7 

(A-1) 

where 

H\v  =  1^;Il2(2R)  +  IIV,'Il2(J?)  • 

Hence,  (A.l)  shows  that  all  elements  ip  of  V  are  bounded  on  JR.  The  limit 

lim  =  0  (A. 2) 

l2~~¥~00 

must,  consequentially,  hold  for  every  (without  loss  of  generality)  positive  real 
number  £  and  every  function  ip  in  V.  Evidently,  then,  such  a  function  ip  can  be 
expressed  as  follows: 

*«.)  =  /  ^  (0(Oe_i(£2^)2)  de  Vx2  6  JR.  (A.3) 

—  OO 

Thus,  from  (A.3)  and  the  Cauchy-Schwarz  inequality  it  is  clear  that 

+oo  +oo 

«*2)|<  /  ^r^W0|e-»(£2^,!d£+  j  |*'({)| 

—  OO  — OO 

—  7M  ^  lr^,|L*(iR))  V(x2,t)  €  JR  X  JR+.  (A.4) 

It  is  then  obvious  from  (A.4)  that 

II^Il~(K)  ^  7!"‘  ^  ^  €  -®+-  (A.5) 

7  See  Aubin  [9]  for  a  demonstration  of  this  fact. 
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Now,  minimize  (A. 5)  with  respect  to  £  to  obtain  (3.4.8).  Note  that  the  constant 
2(-f  )<  in  (3.4.8)  may  not  be  the  sharpest  possible  one  for  an  estimate  of  this  type. 
That  is,  there  may  exist  a  function  in  V  more  optimal  than  the  Gaussian  used 
in  (A.3)-(A.5). 

Appendix  B.  In  this  appendix  the  estimates  (3.4.9)  and  (3.4.10)  for 
IMIlj(2R)  ||<^>/ lz,2(2R)  established.  First  consider  (3.4.9).  From  the  in¬ 
tegral  equation  in  (3.2.15),  the  Minkowski  inequality  and  the  Cauchy-Schwarz 
inequality  it  is  clear  that 


^  ^  llz,2(iR2)  IV?llz,2(iR)  +  ^  ll/slll,2(Ji)  •  (B-l) 

With  the  aid  of  the  decomposition  of  fa  provided  in  (3.3.17),  the  bound  (3.3.21), 
and  the  fact  ||^a lz,2(iR)  =  ^  IKIIl^JR)  (B.l)  implies  that 


II  V^ll  Z.2(  JR)  ^ 


A(7T  + 

1  —  A  \K„\ LHJR2) 


(B.2) 


Now,  use  (3.4.8)  and  (3.4.9)  in  (B.2)  to  give  (3.4.9). 

Next  consider  (3.4.10).  Recall  that  in  order  to  obtain  an  estimate  for  the  L 2 
norm  of  <p'  over  M  it  is  sufficient  to  require  that  s  be  an  element  of  X  fl  X,  where 
X  is  given  by  (3.4.17).  Suppose  that  this  is  the  case.  Then  it  is  permissible  to 
differentiate  the  integral  equation  in  (3.2.15)  to  obtain 


+oo 

*>'  +  A  fK.(;0v(()di  =  lf:  on  JR,  (B.3) 

— oo 


where 

I 

=  2K](z,,()3(Z2)  ~  “p  +  2 L.(Xl,()  V(I2,«6JJ2,  (B.4) 

,  X2  -  ; 


L 
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and 

T  ,  *v  _  s(x2)  -  s( 0  -  (aa  -  CK(g2)  +  i(f2  -  QV(f 2) 

”  (xa  -  0((-(*2)  -  *(0)2  +  (xa  -  O2) 

V(x2,0€iR2  (B.5) 

Clearly,  with  the  aid  of  the  Cauchy-Schwarz  and  Minkowski  inequalities,  (B.3) 
implies  the  following  estimate  for  Yp'Xlur)' 

IV’I  mm.)  ^  ^\K»\L*(R't)  H^Ilj(2R)  +  ^  If  a  II  £,2(ilR)  •  (B-6) 

So,  given  bounds  for  \Ka\,Li{lR'1)  and  ||/j|Il2(2r)  *n  terms  of  the  L2  norms,  over  M , 
of  the  first  three  derivatives  of  s,  (B.6)  will  provide,  in  conjunction  with  (3.4.9),  an 
estimate  for  fly^H A  bound  for  the  L2  norm  of  L„  over  JR2  can  be  obtained 
in  exactly  the  same  manner  as  that  established  for  Ka  in  Section  3.3.  This  bound 
is 

^  ls"l£*(jt)  \s'"\h(]R)  •  (®*^) 

Now,  from  (B.4),  (B.5),  (B.7),  (3.3.2)  the  Cauchy-Schwarz  inequality  and  the 
Minkowski  inequality  one  obtains  the  following  estimate  for  || ||  : 

^  Is,|l~(.R)  15,,Il»(K)  I,s/ ll£,2(iR)  I5,,|£°°(jr) 

+ #  i*"ii(«)  ■  (b.8) 

It  is  also  clear  that 


ll/«  llz,*(]R)  ^  n  I51l*(JR) 

+  I  Is,,|l=o(jJ)  I5,|l*(JR)  •  (B.9) 

Using  (3.4.8),  inequalities  (B.8)  and  (B.9)  become 

^  !-s# IZ*(jr)  I5  II5  l£~(iR) 

+  ^  8-®” tl Z»(jr)  II5”  !£«(«) »  (B.10) 
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and 

I/jIl2(K)  —  77  I5  li2( JR)  +  2(f)'1  U  |£2(/i)  l-S  Ht2(jR)  iz-2(JR2) 

+  (f )  *  ls,|L2(Ji)  I5  IIl2(R)  IIs  li2(R)  l^» Il2(JR2)  >  (B-11) 

respectively.  Combining  (B.6),  (B.10),  (B.ll)  and  (3.3.15)  leads,  after  a  bit  of 
algebra,  to  the  desired  estimate  (3.4.10). 
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Figure  1:  Graph  of  the  shear  stress  response  function  r. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2h  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

Technical  Report  No.  2 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

California  Institute 
of  Technology 


6c  ADDRESS  (City,  State,  and  ZIP  Code) 

Pasadena,  California  91125 


6b.  OFFICE  SYMBOL 
(If  applicable ) 


1b  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

UNLIMITED 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research 


7b  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Pasadena,  California  91106-3212 


9a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


9c.  ADDRESS  (City,  State,  and  ZIP  Code) 

800  N.  Quincy 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

N00014-90- J-1871 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJfc'CT 

TASK 

ELEMENT  NO. 

NO 

NO 

22217-5000 

WORK  UNIT 
ACCESSION  NO 


On  the  Construction  of  Two-Phase  Equilibria  in  a  Nonelliptic  Hyperelastic  Material 


12  PERSONAL  AUTHOR(S) 

Eliot  Fried 


13a.  TYPE  OF  REPORT 

Technical 


16.  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED 
FROM _  TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  115.  PAGE  COUNT 

September,  1990 


[•Tl 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  work  focuses  on  the  construction  of  equilibrated  two-phase  antiplane  shear  deformations 
of  a  non-ell iptic  isotropic  and  incompressible  hyperelastic  material.  It  is  shown  that  this 


_ _ _r_  _  __  - - ’  ~  r  ■  -  -  .  ^  .  w  i  ijr  |7l-i  u  iuj  y  iu  mu  lci  iui.  *  L  ■  O  cmuwn  LIICI  L  Lilli 

material  can  sustain  metastabl e  two-phase  equilibria  which  are  neither  piecewise  homogeneous 
nor  axisymmetric,  but,  rather,  involve  non-planar  interfaces  which  completely  segregate 
inhomogeneously  deformed  material  in  distinct  elliptic  phases.  These  results  are  obtained 
by  studying  a  constrained  boundary  value  problem  involving  an  interface  across  which  the 
deformation  gradient  jumps.  The  boundary  value  problem  is  recast  as  an  integral  equation 
and  conditions  on  the  interface  sufficient  to  guarantee  the  existence  of  a  solution  to  this 
equation  are  obtained.  The  constraints,  which  enforce  the  segregation  of  material  in  the  twi 
elliptic  phases  are  then  studied.  Sufficient  conditions  for  their  satisfaction  are  also 
secured.  These  involve  additional  restrictions  on  the  interface  across  which  the  deformatioi 
gradient  jumps — which,  with  all  restrictions  satisfied,  constitutes  a  phase  boundary.  An 
uncountably  infinite  number  of  such  phase  boundaries  are  shown  to  exists  It  is  demonstrate! 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
0UNCLASSIFIED/UNLIMITEO  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

J.  K.  Knowles 


DO  Form  1473,  JUN  86  i 


21  ABSTRACT  SECURITY  CLASSIFICATION 

□  dtic  users  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(818)  356-4135 


Previous  editions  are  obsolete. 


security  classification  OF  THIS  PAGE 
UNCLASSIFIED 


